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Abstract Long-lived, stable jets are observed in a wide variety of systems, from proto- 
stars, through Galactic compact objects to active galactic nuclei (AGN). Magnetic fields 
play a central role in launching, accelerating, and coUimating the jets through various me- 
dia. The termination of jets in molecular clouds or the interstellar medium deposits enor- 
mous amounts of mechanical energy and momentum, and their interactions with the external 
medium, as well, in many cases, as the radiation processes by which they are observed, are 
intimately connected with the magnetic fields they carry. This review focuses on the prop- 
erties and structures of magnetic fields in long-lived jets, from their launch from rotating 
magnetized young stars, black holes, and their accretion discs, to termination and beyond. 
We compare the results of theory, numerical simulations, and observations of these diverse 
systems and address similarities and differences between relativistic and non-relativistic jets 
in protostellar versus AGN systems. On the observational side, we focus primarily on jets 
driven by AGN because of the strong observational constraints on their magnetic field prop- 
erties, and we discuss the links between the physics of these jets on all scales. 

Keywords magnetic fields ■ jets • polarization ■ radio continuum: galaxies ■ X-rays: 
galaxies 



1 Introduction 



Astrophysical jets have been discovered in a very wide range of physical systems, from re- 
gions of star formation, to gamma-ray bursters and stellar black holes, on up in scale and 
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energy to the nuclei of giant radio galaxies and quasars. They are often associated with ac- 
cretion discs and their launch, internal dynamics, and impact upon their surroundings are 
of ongoing importance and significance over a whole range of astrophysical phenomena 
— from the parsec (pc) scales that characterize the impact of jets from young stellar ob- 
jects within their natal molecular clouds, to the megaparsec (Mpc) scales of giant clusters 
of galaxies whose intercluster media a re significantly stirred a nd heated by jets emanating 
from the cores of active galaxies (e.g.. lDurm and Fabianll2006i) . On stellar scales, jets from 
young stellar objects (YSOs) span the breadth of the stellar mass spectrum: from massive 
stars such as W33A whic h shows a hot fast jet when observed in Br 7 lines as well as an 
underlying accretion disc dPavie s et al.'2010^■ to the well known solar-mass systems such as 
classical T-T auri stars (TTS) (see Ray et al. 2007, for a review), and on down into the brown 
dwarfs (e.g., IWhelan et al.ll200^ . There is good evidence that all of the protostellar sys- 
tems are associated with discs. Stellar mass black holes, otherwise known as microquasars 
(§t/Iirabel 2005) are also very active sources of relativistic jets. On extragalactic scales, jets 
are observed to be associated with accretion on to supermassive bla ck holes (up to a few bil- 
lion solar masses) in the centres of active galaxies (lBlandfordll200lh . Models for gamma-ray 
bursts (GRBs), which are also located at cosmological distances, focus on either the col- 
lapse of massive stars or the merger of two compact stellar objects to form a black hole in 
a process that releases about a solar mass wor th of gravitation al potential energy and which 
generates a highly relativistic magnetized jet (lMeszarosll2012h . 

Given the vast range in energies, scales, and details spanned by these diverse systems, 
it is intriguing that accreting systems seen in nature make such extensive and ubiquitous 
use of energetic outflows as part of their evolution. Is there a universal mechanism at work 
that almost guarantees that jets will be present? What is the basic physical mechanism that 
creates and drives jets across all of these scales? What role do they play in the evolution 
of such accreting compact systems? The fact that all these systems share common physical 
components — central gravitating sources, accretion discs, and magnetic fields — provides 
the unifying foundation on which theoretical models are constructed and observationally 
tested. 

Magnetic fields play a critical role in creating, driving, and directing the evolution of jets. 
Indeed, the successful existing models of jets in these diverse systems invoke magnetic fields 
threading central objects such as gaseous accretion discs or rapidly rotating central objects 
themselves (see reviews by, e.g., , Shu et al. 2000; Pudritz et aL 2007). The release of gravi- 
tational potential energy during accretion is a key power source for most of these jets. The 
great utility of magnetic torques is that they can readily tap into the shear energy available 
in the accretion discs that form as a consequence of gravitational collapse. Models of mag - 
netized outflows from accretion discs in AGN, developed by iBlandford and Pavnd 1 1982h. 
were applie d soon after to understanding jets in protostellar svstems IPudritz and Normmil 
Il98lll986h . These theories envisage that material from a disc, in being forced to co-rotate 
with a threading magnetic field, is flung outwards along sufficiently inclined, open field 
lines. The gas continues to accelerate as long as the field is strong enough to maintain co- 
rotation with the central disc, or rotor. The inertia of this matter moving along the field 
line causes the field to wrap up in a toroidal-like structure — creating a helical field. It 
is this dominant, outer toroidal field component that exerts a pinching force upon the jet, 
collimating it towards the axis. A second major prediction of this basic theory is that jets 
should be extracting the angular momentum of their sources, and jet rotation is a key pre- 
diction of the models. Another abundant source of energy resides in the rotation available 
in rapidly rotating black holes or pulsars, a nd once again it is known th at electromagnetic 
torques can tap rotational energy efficiently dBlandford and Znaiekll 1977b . So for rapidly ac- 
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creting or spinning objects, large scale magnetic fields appear to be nature's preferred tool 
in extracting energy and angular momentum and directing them in the formation of ordered, 
magnetically coUimated, outflows. 

Observations of crucial properties of jets, such as their densities, temperatures, velocity 
structure, and magnetic fields, vary quite considerably between protostellar and extragalac- 
tic systems. In the former, forbidden line emission (FLER) produced by shocks in the jet 
allows one to measure directly everything with the exception of the magnetic field. On the 
other hand, extragalactic jets emit synchrotron radiation which arises from the relativistic 
electrons that spiral around magnetic field lines. This radiation allows one directly to mea- 
sure the direction of the magnetic field in the jet, and to estimate its strength. However, the 
absence of distinct emission lines from the outflowing material itself in AGN jets prevents 
one from measuring the other crucial physical variables. It is for this reason that protostellar 
systems have become such important laboratories for jet studies. Until 2010, these limita- 
tions prevented direct observational comparison of these classes of jets, and, as a result, 
meant that it was impossible to test whether or not jets have a universal magnetic field based 
mechanism. This all c hanged with the discovery of s ynchrotron radiation from the proto- 
stellar jet, HH 80-81 Ic arrasco-Gonzalez et al.ll2O10l) . Relativistic electrons in otherwise 
non-relativistic jets arise from acceleration in the jet's terminal shock. The results show that 
the magnetic structures of protostellar and AGN jets are similar — with a field parallel to 
the jet axis at the centre and a surrounding helical field towards the jet edges. 

In this paper, we examine the structure and role of magnetic fields in the physics of jets 
— from their launch on discs or compact rotors, to their dynamics and evolution, and finally 
to a consideration of how such strongly magnetized flows affect their environments on these 
various scales. We begin our analysis with non-relativistic systems — the protostellar jets 
that are a ubiquitous and vital component of any star-forming region (Section 2). We also 
summarize the theory and simulations of relativistic jets in the context of AGNs. (We do 
not consider the impulsive, short-lived jets in gamma-ray bursts, although much of the same 
physics is expected to ap ply, as the ob serva tional challeng es are very different; the reader 



relativistic systems and discuss the measurement of the structure and dynamics of AGN jets 
on highly compact (pc) scales, probing the conditions near at the source of the jets using 
VLBI. The magnetic structure of jets on these scales is probably a reflection of their internal 
dynamics and not of the shocks which occur throughout the body of the jets. Indeed, the 
mounting evidence for helical magnetic fields in extragalactic jets, taken with the newly 
discovered evidence for helical fields in protostellar jets discussed above, strongly suggest 
that the physical mechanism of coUimation may be the same. Finally, we go out to the 
scales of jet termination (Section 4), focussing on the evidence for magnetic field properties 
on these largest scales and exploring the connection between magnetic fields in jets in the 
hearts of AGNs with jets at the largest scales (1 kpc-1 Mpc). We conclude our study with a 
discussion of future prospects in this field (Section 5). 

2 Magnetic launching of jets 

2. 1 Observed properties of protostellar jets 

Observations of the launch regions even of protostellar jets face significant observational 
challenges. For the nearest regions of star formation and jets in the Taurus molecular cloud 
(distance 140 pc), l" resolution corresponds to a physical scale of 140 AU. Theoretical 




Section 3 we then turn to 
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models discussed below picture a jet acceleration region which is of the order of several tens 
of AU (O.l"), and the launch region in the very inner regions of the disc (several AU) that 
may be an order of magnitude smaller (0. 01"). Jets around young stars have been studied 
on scales ranging from tens of AU to 5 pc jBallv et al]l2007h . 

Hubble Space Telescope {HST) observations of jets from YSOs have advanced our ex- 
perimental understanding of jets enormously over the last decade, probing physical scales 
down to tens of AU. There are four basic types of observations that point towards the strong 
coupling of jets with underlying accretion discs — the correlation of jets with discs, the 
onion-like velocity structure of jets, jet rotation, and the strong link between the rates of 
mass transport in the jet to the accretion rate in the underlying disc. The advent of high res- 
olution, spectro-imaging of jets by the HST and ground-based observatories equipped with 
adaptive optics has revealed their remarkable internal structure. Jets are strongly correlated 
with the presence of protostellar accretion discs around their host stars — at least during the 
first million years or so of their pre-main sequence evolution. Protostellar jet properties can 
be directly probed by observations of the forbidden line emission that they produce. Optical 
forbidden lines of oxygen, sulphur and nitrogen, for ex ample, are produced in very specific 
ranges of density and temperature (see lRav et al.ll2007l for a review). The spectro-imaging 
results indicate that jets have an "onion-like" velocity structure, with the highest-speed out- 
flow component, up to 400 km s^^ in some sources, found at the core of the jet, while 
the lower-velocity components are found systematically farther from the outflow axis. The 
highest-resolution (HST) studies of the launch region of TTS jets (about 50 AU) show gas 
densities > 2 x 10'' c m""^, modest electro n temperatures (T ~ 2 x 10"* K), and low ioniza- 
tion levels (0.03 - 0.3) dCoffev et al.l2008h . As we shall see, the onion-like velocity structure 
seen close to the source is a reflection of the differential rotation of the underlying Keple- 
rian disc. Higher-resolution observations, down to the 1-AU scales predicted for jet launch, 
would be highly desirable. 

The discovery of jet rotation is perhaps one the most significant advances of the last 
decade jBacciotti et aLll2002l : ICoffev et alj|2004l) . By orienting the slit of the spectrograph 
perpendicular to the jet axis at some distance along the jet, velocity asymmetry that is com- 
patible with jet rotation can be observed. Shifts range from 5-25 km s^' at positions of 25 
AU along the jet. From the point of view of jet dynamics, these observations allow one to 
measure the amount of angular momentum that is being carried by the jet — which ranges 
from ~ 60% to ~ 100% of the angular momentum that would be carried through the disc by 
the observed disc accretion rates. The ratio of mass carri ed in the jet to that moving through 
the underlying accretion disc, Mj^t/Ma 0.01 - 0.07 dCoffev et al.1l2008h . 

The strength of the magnetic field in the HH80-81 jet is estimated from the synchrotron 
emission (on the basis of equip artition arguments: see Section 14.2b , to be 0.2 milligauss 
jCarrasco-Gonzalez et al.ll2010l) . Even though the temperatures of jets are quite low, the 
mechanism of diffusive shock acceleration has been shown to be able to accelerate electrons 
moving at a few hundred km s~' up to relativistic speeds. If this technique can be exploited, 
we may indeed be on the threshold of being able to measure all of the basic physical quan- 
tities that are needed to test the physics of protostellar jets. 

2.2 Structure and dynamics of non-relativistic jets 

At the heart of the question of why jets are ubiquitous is the fact that rapidly rotating mag- 
netized bodies undergo strong spin-down torques exerted by the threading fields. The basic 
theory of magnetized winds from rapidly rotating stars was first worked out in the seminal 
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paper bv lMestell ( fl96^ . This theory focuses on stationary solutions of axisymmetric flows 
from magnetized bodies — as applied in particular to rotating stars. It has turned out that this 
theoretical framework is very useful in analyzing the structure of jets and of the structure of 
their magnetic fields. 

The overall dynamics of jets, and the physics of how they are launched, is complicated. 
However, significant progress has been made by the application of the theory of stationary, 
2D (axisymmetric) MHD flows. The usefulness of this highly simplified treatment arises 
from the fact that we can derive conservation laws that turn out to be important even for un- 
derstanding the behaviour of stationary jets in 3D. First we decompose vector quantities into 
poloidal and toroidal components (e.g. magnetic field B = Bp + B^e^). In axisymmetric 
conditions, the poloidal field Bp can be derived from a single scalar potential a(r, z) whose 
individual values, a = const, define the surfa ces of constant magnetic fl ux in the outflow 
and can be specified at the surface of the disc (Pelletier and Pudritall992h . 

The joint conservation of mass and magnetic flux along a field line can be combined 
into a single function k that is called the "mass load" of the wind, which is a constant along 
a magnetic field line: 

pvp = fcBp. (1) 

This function represents the mass load per unit time, per unit magnetic flux of the wind. For 
axisymmetric flows, its value is preserved on each ring of field lines emanating from the 
accretion disc. Its value on each field line is determined by physical conditions — including 
dissipative processes — near the surface of the rotor. 

The mass load plays a central role in jet dynamics (see Sections 3 and 4) and controls jet 
rotation, coUimation, and angular' momentum extraction. It may be more revealingly written 
as 

%) = a = ^^, 

where dMw is the mass flow rate through an annulus of cross-sectional area dA through the 
wind and dW is the amount of poloidal magnetic flux threading through this same annulus. 
The mass load profile (as a function of the footpoint radius ro of the wind on the disc) 
is determined by the physics of the underlying disc and imposes an important boundary 
condition for all aspects of jet physics. 

The conservation of angular momentum along each field line leads to the conserved 
angular momentum per unit mass; 

1(a) = rvj, — —7- = const. (3) 

^ ' ^ ink 

The form for / reveals that the total angular momentum is carried by both the rotating gas 
(first term) as well by the twisted field (second term), the relative proportion being deter- 
mined by the mass load. 

The value of l{a) that is transported along each field line is fixed by the position of the 
Alfven point in the flow, where the poloidal flow speed reaches the Alfven speed for the first 
time (mA = 1). Hence, the value of the specific angular momentum is found to be 

l{a) = norl = {TA/roflo. (4) 



where Iq = VKfifo = ^^o^'o is the specific angular momentum of a Keplerian disc. (Note, 
however, that this treatment can easily be generalized to sub-Keplerian discs.) For a field 
line starting at a point ro on the rotor (disc in our case), the Alfven radius is rA('"o) and 
this constitutes a lever arm for the flow. The result shows that the angular momentum per 
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unit mass that is being extracted from the disc by the outflow is a factor of {rj^/ra)'^ greater 
than it is for gas in the disc. For typical lever arms, one particle in the outflow can carry the 
angular momentum of ten of its fellows left behind in the disc. 

The conservation of energy along a field line is expressed as a generalized version of 
Bernoulli's theorem. Since the terminal speed vp = Voo of the wind is much greater than its 
rotational speed, and for cold flows, the pressure may also be ignored, one finds the result: 

Woo ^ 2^/^J7o7'A = {rA/rQ)vescfl- (5) 

There are several important consequences for jet dynamics that follow from this simple 
scaling. The first is that the terminal speed exceeds the local escape speed from its launch 
point on an accretion disc by the lever arm ratio. Furthermore, the scaling predicts that the 
terminal speed scales with the Kepler speed as a function of radius, so that the flow will 
have an onion-like layering of velocities, with the largest velocities inside and the smallest 
on larger scales. As we have already noted, such velocity structure is actually observed in 
the dynamical structure of well-resolved jets. Finally, the terminal speed depends on the 
depth of the local gravitational well at the footpoint of the flow (the dependence on the local 
Keplerian rotation speed in the case of a disc rotor) — implying that it is essentially scalable 
to flows from discs around YSOs of any mass and therefore universal. 



2.3 Determining the jet launch region 

A very useful combination of the energy and angular momentum con servation that we de- 
rived above is given by the Jacobi constant along a field line (e.g., Pelletier and Pu drit3 
[l992): j{a) = e(a) — nol{a). This expression has been used (Anderson et al. 2003) to 
constrain the launch region on the underlying Keplerian disc. Observers measure the jet 
rotation speed, u^^oo at a radius Voo, and the poloidal speed of the jet at this point, up,oo. 
Evaluating j for a cold jet at infinity and noting that its value (calculated at the foot point) 
is j{ao) = — (3/2)i'^ Q, one can solve for the Kepler rotation at the point on the disc where 
this flow was launched: 

(6) 

When this relation is appl ied to the observed r otation of the Large Velocity Component 
(LVC) of the jet DG Tau ( iBacciotti et alJl2002l) . this yields a range of disc radii for the 
observed rotating material in the range of disc radii, 0.3^ AU, and the magnetic lever arm 
is rx/ro — 1.8-2.6. These results show that an extended region of the disc can give rise to 
the jet, and that, moreover, the lever arm is in the predicted range for the efficient launching 
of jets. 



2.4 The disc/jet connection: magnetic torques on discs 

The magnetized wind that is accelerated off the rotor carries away the rotor's angular mo- 
mentum. In the case of an accretion disc, the angular momentum equation for the accretion 
disc undergoing a purely external magnetic torque (i.e., viscous torque neglected) may be 
written: 

• d(rouo) 2 n n I 

dn) ^ "''O-S^-S^lr-o.ff. (7) 

where we have ignored transport by MRI turbulence or spiral waves. Note that the toroidal 
field plays a fundamental role in disc structure and accretion because it is so central to the 
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action of a magnetic torque upon the disc. By using the relation between the poloidal field 
and outflow on the one hand, as well as the link between the toroidal field and rotation of 
the disc on the other, the angular momentum equation for the disc yields one of the most 
profound scaling relations in disc wind theory, namely, the link between disc accretion and 
mass outflow rate: 

M^~{rA/rafM^. (8) 

The observationally well-known result that, in many systems, jMa, ~ 0.1 is a conse- 
quence of the fact that lever arms are often found in numerical and theoretical work to be 
rA/^o ^ 3 - the observations of DG Tau being a perfect example. Finally, we note that the 
angular momentum that is observed to be carried by these rotating flows (e.g ., in DG Tau) is 
a con sistent fraction of the excess disc angular momentum - from 60-100% feacciotti et al.l 
|2004) . which is consistent with the high extraction efficiency discussed here. 

In reality, the evolution of discs is dictated by both viscous as well as wind torques. 
The magnetic torque exerted on the underlying disc has important effects upon its structure 
( Ferre ira 2002; Ferreira and Casse 2004). Compared to standard a-disk models. Combet and Ferreiral 
1 20081) show that outflow torques modify the underlying radial structure of discs, resulting 
in column density jumps of a couple of orders of magnitude in going from the outer SAD 
solution to the inner-jet dominated zone of the disc. The jet-dominated region is also cooler 
and thirmer. 



2.5 Jet collimation and toroidal fields 

In the standard picture of hydromagnetic winds, collimation of an outflow occurs because 
of the increasing toroidal magnetic field in the flow resulting from the inertia of the gas. 
Beyond the Alfven surface, the ratio of the toroidal field to the poloidal field in the jet 
is of the order B^/ Bp ~ r /r^ 2> 1, so that the field becomes highly toroidal beyond the 
Alfven surface. Collimation is achieved by the tension force associated with the toroidal field 
which leads to a radially inwards directed component of the Lorentz force (or "z-pinch"); 

^Lorentz.r — jzB^p. 

The current carried by a jet can be easily written down: 

I = 2tv f j,{r', z')dr' = {c/2)rB^. (9) 
Jo 

The link between magnetic field structure and the collimation of jets was made bv lHewaerts and NormarJ 
11989), where it was shown that two types of solution are possible, depending upon the 
asymptotic behaviour of the total current intensity in the jet. In the limit that / — > as 
r ^ oo, the field lines are paraboloids which fill space. On the other hand, if the current is 
finite in this limit, then the flow is coUimated to cylinders. The collimation of a jet therefore 
depends upon its current distribution — and hence on the radial distribution of its toroidal 

field. 

It can be shown jPudritz et al1l2006l) that, for a power-law distribution of the magnetic 
field in the disc, B2(ro,0) oc ^q"^, with an injection speed at the base of a (polytropic) 
corona that scales as the Kepler speed, that the mass load takes the form k oc r^^^~'"'^\ In 

this regime, the current takes the form I{r, z) oc t-q Thus, the current goes to zero 

for models with ^ < — 1/2, and these therefore must be wide angle flows. For models with 
/i > — 1/2, however, the current diverges, and the flow should coUimate to cylinders. 
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These results predict that jets should show different degrees of coUimation, depending 
on how they are mass loaded. As an example, neither the highly centrally concentrated, 
magnetic field lines associated with the initial split-monopole magnetic configuration used 
in simulatio ns bvlRomanova et alj ( il997D , nor the similar field structure invoked in the X- 
wind model dShu et alJbOOol) ^ should become coUimated in this picture. On the other hand, 
less centr ally (radially) concentrated mag netic configurations such as the potential config- 
uration of lOuved and Pudritj (11997*) and iBlandford and Pavnel ( Il982h should coUimate to 
cylinders. By varying both the density profile of the disc as well as its magnetic structure, 
iFendtl (l2006h showed simulations in which highly coUimated outflows from a "flat" radial 
disk structure produced unsteady knots in the jet. 

The magnetic coUimation picture, although it h as considerable support from theory and 
simulations, has been questioned bv lSpmitlbOlOl) . who notes that the toroidal magnetic 
pressure in jets will force them to expand. CoUimation in this picture results from magnetic 
fields and pressure in the external medium. Simulations do indeed provide evidence for a 
regime in which magnetic bubbles associated with simulations of protostellar outflows are 
observed to expand rad ially. These cases a rise in simulations of disc winds that produce 
weak toroidal jet fields jSeifried et al.ll2012l) . Strong toroidal fields in these jets give rise to 
well coUimated flows. 



2.6 Computer simulations: magnetic field structure of jets and their stability 

Numerical simulations have played an enormous role in exploring the origin and evolution 
of jets. The pioneering stage of simulations of jets from discs took place in the mid 1980s 
through the 1990s, with the advent of magnetohydrodynamics (MHD) codes such as ZEUS. 
These early simulations recognized that, since the mass load of the disc outflow is a key 
variable that controls jet dynamics and structure, then the disc can be taken as a bound- 
ary condition for the flow, with the initial magnetic field structure and mass loading defined 
across this boundary. Simulations demonstrated most aspects of the theory outlined above — 
the coUimation of the jet, the concentration of jet density towards the axis of the flow, the ac- 
celeration to high speeds, the "onion-like" velocity st ructure of the jet, and the dominance of 
the toroidal magnetic field structure can b e seen (e .g.,'Ouved et al.'l 997l : l0uved and Pudritj 
[l997; Romanova et al, 1997; Krasn opolskv et alJll999 : Fendt 200^Highly e pisodic be- 
havio ur of jets arises as a consequence of sufficiently low mass loading of the jets jOuved et al.l 
Il997h . 

Figure [l] show s ZEUS 3D simulations of a jet being launched from an underlying disc 
dStaff et alj|2010h . The two panels show snapshots of the magne tic field line geome t ry and 
jet density that arises with i n two particular jet models — those of lOuved and Pudrita ( Il997h 
and IBlandford and Pavnel ( Il982h . The field lines have a distinct toroidal component that 
dominates the structure of the overall field on larger scales. Poloidal magnetic field dom- 
inates towards the axis of the outflow. The bow shock created by the jet is clearly seen in 
both cases. The general structure of these models is rather similar. The densest part of the 
jet are surrounded by a diffuse cavity-like region that is dominated by the toroidal magnetic 
field. The more extended outer material moves more slowly than the jet core. 

If jets are dominated by toroidal fields far from the jet axis, why are they stable? The- 
oretical studies of the linear stability of jets show that they should become unstable when 
their speeds exceed the local Alfven velocity (Ray 1981). This was further demonstrated by 
3D simulations of idealized radio jets by iHardee and RosenI ( Il999i) . It turns out, however, 
that toroidally dominated jets stabilize themselves through nonlinear processes. Simulations 
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Fig. 1 Density and magnetic field lines for the Ouyed and Pudritz (1997) and Blandford and Payne (1982) 
models for the distribution of disc magnetic fields. The snapshots are shown as the jets have propagated out to 
60 AU in the medium surrounding the magnetized disc. The protostar and the accretion disc are in the lower 
left comers of the panels, and hidden by the dense disc corona. Note that the core of the jet is dominat ed by 
poloidal field lines, while the outer regions are predominantly toroidal. Adapted from lStaff et al.|j2010l) . 



show that jet instability, particularly in the form of m = 1 helical modes, builds up beyond 
the Alfven surface lOuved et al.. .2003i) . These are stabilized by the magnetic tension that 
stretches field in the linearly unstable region, which ultimately increases the local Alfven 
speed again. The net result is a jet that is non-linearly stabilized by this mechanism, as well 
as by the "backbone" provided by the poloidal field that dominates near the outflow axis. 

Since jets and accretion discs are highly coupled, jet torque and evolution will directly 
affect the disc. Hence, there has been a major effort over t he last decade to simulate jets and 
their discs as a single dynamical syste m (e.g.. lKudoh et a l. 2002; Casse and Keppens 2 0o3 : 
Ivon Rekowski and Brandenburdl2004h . The important advantage of this approach is that one 
can avoid ad hoc descriptions of the mass loading and disc magnetic field line structure since 
these are self-consistently computed. These simulations show that jet structure and dynam- 
ics retain th e basic features that the simpler simulations capture. Moreover, simulations by 
IZarmi et aTl C2007 ) which include the disk, and follow the evolution of jets using the Adap- 
tive Mesh Refinement code FLASH, find that more than 90% of the gravitational potential 
energy liberated in the accretion flow is released into the jet. Thus jets are highly efficient in 
extracting both energy and angular momentum from the underlying accreting systems. 

2.7 Jets from magnetized star-disc systems 

So far we have emphasized disks as the underlying magnetized engines for jets. However, 
coUimated outflows may derive from both the central stars as well as their surrounding d isks. 
T-Tauri stars (TTS) are well known to be significantly magnetized djohns-KrulllbOOTh . As 
evidence that this may be occurring, we focus on the fact that TTS are also known to rotate 
very slowly for the amount of angular momentum that they are known to be accreting. This 
has long been known to imply that there must be some external torque being exerted upon 
the rotating star, canying off or intercepting the angular momentum. It has been proposed 
that this torque takes the form of an interaction of the star's magnetosphere and the sur- 
rounding disk through closed magnetic field lines (lKonigllll991l) . Recent work suggests that 
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most of the field in such interaction is readily opened up by shear i n the disk, leading to ver y 
little spin-down torque and connection between star and disk (e.g.. lMatt and Pudritzl (l2004h . 
The interaction between inner disk and stellar magnetosphere has also been hypothesized to 
give rise to a so-called "X-wind", wherein the origin of jets and outflows is pictured to be 
an MHD outflow — based o n the principles w e discussed above — and originating from the 
very inner edge of the disk dShu et alj|2000h . In this picture, dipolar stellar magnetic field 
lines interacting with the inner edge of the disc are opened, leading to an outflow which 
intercepts angular momentum flowing inwards through the disc and driving it out of the sys- 
tem in the outflow. One difficulty with this models is in accounting for the huge amount of 
angular momentum that jets are observed to transport — 60% or more of the angular mo- 
mentum transport through the disc at scales of 2 AU or so (Coffey et al. 2008); another issue 
is the difference be tween observed ma gnetic field strengths and the field strengths needed to 
drive the X-winds (|johns-Krul]||2007l) . A third class of models has been proposed, wherein 
an accretion-po wered wind from a magnetized, central star carries off the disk angular mo- 
mentum dMatt and Pudritzll2005b . The energy driving this wind is ultimately drawn from the 
gravitational energy that is released by the infalling material from the disc as it falls along 
magnetospheric field lines connecting the star and the disc. A strong flux of Alfven waves 
stimulated b y the impact of the accretion flow upon the star creates a wave-driven wind 
dMatt and Pudritz 2005 . 2008 ; Cranmer 2008) . 

Simulations of such combined magnetized star-disk systems have been carried out in 
which both the central rotating magnetized star and the surrounding disk are treated as 
boundary conditions for the flows (Fendt 2009). Compared to the case of pure disk winds, 
the overall outflow is initially less collimated. This changes with time, however, as a highly 
coUimated flow emerges as a new dynamical state across the grid when a stellar wind is 
included. The magnetic interaction between the two wind components gives rise to recon- 
nection and large-scale flares as well, and the time-scale for the flaring is long — of the order 
of hundreds of inner disc rotation periods. Finally, if the outflow is dominated by the stellar 
wind, the outflow is either too weak or too high for low and high mass loads respectively. 
The disk jet is required to have stable collimated flow out to large scales. 

The interaction of spinning, magnetized stars wit h their surrounding disks have been ex- 
tensiv ely simulated for axisymmetric configurations dRomanova et alJ 20021 ; iBessolaz et alJ 
l2008h . Magnetospheric accretion has been likened to "funnel flow" wherein matter falls 
onto the star along field lines that are connected to the star and co-rotate with it. The matter 
moves due to the gravitational force, which dominates over the centrifugal force. The lift- 
ing of material off of the disc and the onset of the funnel flow starts at the radius at which 
the disk is slowed due to interaction with the more slowly spinning magnetosphere. These 
simulations showed that the interaction is time-dependent, with bursts of accretion followed 
by quiescent periods in which material from the disc piles up against the magnetosphere. 
Models usually assume that the disks are characterized by a model "a-disc" viscosity. In 
self-consistent 3D simulations, accretion flow from a turbulent disk onto a tilted magnetized 
star, w here turbulence is driven by the MRI instability, has been modelled dRomanova et alj 
l2012h . The pattern of alternating compression due to accretion and reconnection is also seen 
in these 3D simulations. The magnetosphere truncates the disk at a few stellar radii. The 
magnetic structure of the disk outside of the stellar magnetosphere has a toroidal field 
that is a few times larger than the poloidal field component. As we have noted several times, 
the transport of angular momentum by twisted field lines is an important part of the outflow 
mechanism. In the star-disc interaction, the magnetic stresses dominate the matter stresses in 
the disk and the pattern is quite inhomogeneus. This goes together well with the predictions 
outlined in the basic theory above. Finally, recent simulations of star-disc systems report the 
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launch of fast outflow at the interface in systems undergoing very high accr etion rates — 
wherein the magnetic field in the star is strongly compressed dLii et al.ll2012h . No X-winds 
are observed for the much lower accretion rates that dominate most of the star's accretion 
teessolaz et al]|2008l) . 



2.8 Discs and jet formation during gravitational collapse 

We now leave the domain of stationary outflows, and examine how they arise in the fully 
dynamical environments that characterize the birth of stars. Protostellar jets appear during 
the earliest stages of star formation and are the first obvious manifestation that gravitational 
collapse and star formation are underway. No 3D simulations of purely hydrodynamic col- 
lapse show that sustained jets can be produced in collapsing systems. This may be why 
they were missed in early hydrodynamical simulations of star formation. Add a threading 
magnetic field, however, and the launch of energetic jets and outflows is readily observed. 
The formation and evolution of jets takes place in a few related phases: (i) the formation 
and early evolution of a rotating, magnetized dense core, which is usually found embedded 
within larger scale filaments within molecular clouds; (ii) the extraction of angular momen- 
tum from the rotating core by torsional Alfven waves during the early formation and collapse 
stages, which significantly de-spins the core; (iii) the gravitational collapse of a magnetized 
core with the ensuing formation of a protostellar accretion disc; (iv) the creation and launch 
of an outflow from the disc as it forms; and (v) the accretion of the disc onto the central star 
and with it, the gradual disappearance of the jet. 

An important measure of the magnetization of the magnetized core is given by the so- 
called mass-to-flux ratio, which is the ratio of the gravitational to the magnetic energy of 
the system. Systems with supercritical ratios (values greater than unity ) cannot be supporte d 



u pporte d 



against collapse by their threading fields. The ratio can be written as (iMouschovia; 

crit 

A second parameter of importance is the ratio of the rotational to the thermal energy of the 
core, known as /3rot. The phase space spanned by these two parameters turns out to divide 
up the physical properties of discs. 

During the gravitational collapse of magnetized cores, the magnetic field threading the 
region is wrapped up — much like a twisting rubber band. As this happens, an outward 
flow of torsional Alfven waves is created which transports angular momentum. Analytic 
solutions of the braking torque that arises for idealized, disc-like rotors have been written 
down. These have been found to fit the early stages of simu lations of the collapse of a 
magnetized, Bonner-Ebert sphere (iBaneriee and F^idritzlbood) . As the collapse picks up, 
the winding of the field becomes particularly strong at the accretion shock. It is from this 
region that the first strong aspects of outflow begins. Soon after this happens, an irmer region 
of outflow is launched which is centrifugally driven. 

For fields that are sufficiently supercritical, the magnetic braking can be so significant 
that the disc material is left rotating at significantly sub-Keplerian levels. This has raised the 
question of whether or not discs can even f orm in th e reasonably magnetized environm ents 
that a re found in regions of star formation dHermeb elle and Fromang 2008; Mellon an d Lil 
l2008h . The results depend quite sensitively on the initial conditions chosen — i.e., the mass- 
to-flux ratio. At values of ~ 2 - 10 , the resulting systems are strongly affected. This 
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intriguing result points to the importance of having appropriate astrophysical initial condi- 
tions. Recent simulations that feature long-term evolution of discs using "sink" particles find 
that in fact a reasonably Keplerian initial disc can form during the long-term evolution of the 
collapse (Duffin et al., in preparation). A grid of models for the formation of massive stars 
show s that reasonably qui ckly rotating cores and lower magnetized cores produce Keplerian 
discs ( ISeifriedetall201lh . 

Figure 12] shows a zoomed-in picture of the magnetic structure that arises in the collapse 
of a solar mass system of a low mass star formation (Duffin et al. 2012, in preparation). The 
jet produced in this simulation reaches out to the extent of the initial Bonner-Ebert sphere, 
of about 10^ AU in diameter. The gravitational collapse has pulled in the field lines which 
thread the disc. A transient, large disc of 2000 AU in extent has formed, which ultimately 
fragments into a ring and an inner 100 AU scale disc within an outer disc rotation period. 
The irmermost region of the disc is seen to warp, which is a consequence of the torques 
induced by the magnetized wind. The outflow, as a consequence, precesses with the warp 
in the disc. Precessing jets are common and the simulation represented in this image shows, 
for the first time, that magnetic fields can cause this behaviour. 

The launch mechanism of the outflows during collapse is under some current debate. For 
low-mass star formation (A/core ~ IMq), the generation, evolution an d properties of proto- 



stella r outflows have been s t udied in great detail over the last years (e.g., Baneriee and Pudrit3 
l2006l : iMellon and Lill2008l : iHennebelle and Fromand l2008l : louffin and Pudritzll2009l) . The 



two basic mechanisms that have been discussed in the literature focus on whethe r the out- 
flow i s driven by centrifug al acceleration (Blandford and Payne 1982; Pudritz and Nomiaril 
1 19861 : IPelletier and Pudrit3ll992l : lFerreiralll997h or by the pressure of th e toroidal magnetic 
field that was wrapped up in the initial collapse jLvnden-Belll 1 199^ , [20031) . The second 
mechanism — known as a tower flow — is highly transient and does not lead to sustained 
outflow moving out from the source. It is based on the study of equilibria of highly wound 
magnetic structures. General energy theorems demonstrate that they form tall magnetic tow- 
ers, the height of which grows with every turn at a velocity related to the circular velocity in 
the accretion disc. The pinch effect amplifies the magnetic pressures toward the axis of the 
towers. 

A more generalized criterion for outflow in disc wind theory , assuming that the gas is 
not strictly co-rotating with the field, has recently been proposed ( ISeifriedetal.1l2012h . This 
model shows that disc-wind theories can have acceleration mechanisms that are driven by 
centrifugal forces as well as toroidal fields. Both magnetic tower and disc-wind theories can 
be commonly understood as a consequence of the same disc- wind models. It is of particular 
interest that the MHD flow equations discussed above actually have two aspects to them. 
Close to the outflow axis, it may be shown that the toroidal field generation is not significant, 
and there the outflow is driven like a centrifuge. Farther away, however, a toroidal field 
component is inescapably produced. It can be shown that the pressure gradient associated 
with this contributes to pushing a continuous outflow from the outer regions of the disc. 

The two regions of the jet may be seen in Fig.[3]where we compare the size of the region 
controlled by purely centrifugal effects with those arising from toroidal pressure gradients. 
This illustrates that magnetic fields in jets play different roles depending upon where they 
are found. The difference between these two snapshots is the region in the jet where gas is 
driven primarily by toroidal field gradients. Note that these dominate in the outer parts of 
the jet, whereas the centrifugal driving predominates in the region close to the jet axis. 
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Fig. 2 Magnetic field structure in a collapsing magnetized core. A flattened structure is produced by the 
gravitational collapse. The purple contours correspond to 3.3 X 10^^^ g cm^'' (n = 8.5 X 10*^ cm^'') and 
black contours to 1.3 X 10^^^ g cm~'^ (n = 3.4 X 10® cm~'^). Yellow contours coiTespond to outflowing 
velocities Vz > 1.5 km s^^, and extend to a height of 10^ AU. Black hues are magnetic field lines. Adapted 
from Duffinetal. (2012). 



2.9 Relativistic considerations and relation to protostellar jets 

There is nothing in the physics of jets that suggests that relativistic sources, such as micro- 
quasars and AGN, and should have radically different driving mechanisms from those of 
protostellar systems — indeed, the dominant magnetic models make use of similar MHD 
equations. The fact that both a sufficiently rapidly rotating black hole and/or the surrounding 
magnetized disc are sources for relativistic jets has an obvious parallel with the protostel- 
lar systems. The theory and simulations of non-relativistic jets, as applied to protostellar 
systems, show some remarkable connections with the observations of AGN jets, which are 
discussed in Section 3 and Section 4. There are some important differences, however, which 
we discuss here. 
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Fig. 3 Left: the region of jet that is driven by centrifugal acceleration — denoted in grey — in a simulation 
of a collapsing, solar mass, magnetized Bonner-Ebert sphere. Right: as left, but grey regions show the regions 
controlled by the full acceleration mechanism, including toroidal pressure gradients. Adapted from Duffin et 
al. (2012). 



The th eory of relativistic jets shows th at the essential parameter is the magnetization 
parameter jMichelll 19691 : ICamenzindl[l98^ , 

(11) 



4Mc3 



where the iso-rotation parameter Hp is the angular velocity of the magnetic field lines. The 
function <]/ = Bpr^ is a measure of the magnetic field distribution and M = upvpR^ is 
the mass flow rate within the flux surface. Equation (lllb demonstrates that the launch of 
a highly relativistic jet requires at least one of three conditions — rapid rotation, a strong 
magnetic field, and/or a comparatively low mass load. A strong magnetic field lowers the 
mass load in the jet, with the consequence that more Poynting flux can be converted into 
kinetic energy per unit mass flux. 

In the case of a spherical outflow = constant) with negligible gas pressure one may 
derive the Mic hel scaling between the asymptotic Lorentz factor and the flow magnetization 
jMichelll 1969h . 

Too = a^/^ (12) 

Depending on the exact magnetic field distribution tf'(r, z), in a c ollimating jet the matter 
can be substa ntially accelerated beyond the fast magnetosonic point ( Beg elman and Lill994 
iFendt and Camenzind.19 9 6.). As a result, the power law index in eg. ([T2ll can b e different 



from the Michel-scaling jFendt and Camenzindl[l996l : Iviahakis and K6nigj[2003l ) . 



An essential difference between relativistic and non-relativistic flows is the role of 
electiic fields. They are negligible in non-relativistic flows but are comparable to mag- 
netic fields in the relativistic case. The issue is that these fields can lead to the decoUi- 
mation of the jet. Current-carrying relativistic jets have, however, been demonstrated to col- 
limate (Chiueh et al. 1991). Relativistic jets have another feature with no counterpart in non- 
relativistic systems, namely the light cylinder, located at the cylindrical radius rj = c/ fl-p. 
At the light cylinder the velocity of the magnetic field lines "rotating" with angular velocity 
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i7p coincides with the speed of light. Outside the Ught cylinder, the magnetic field lines 
"rotate" faster than the speed of light. As the field line is not a physical object, the laws of 
physics are not violated. The light cylinder is interpreted as the Alfven surface in the limit 
of vanishing matter density (force-free limit). Crucially, the location of the light cylinder de- 
termines the relativistic character of the magneto sphere. If the light cylinder is comparable 
to the dimensions of the object investigated, a relativistic treatment ofMHD is required. 

Simulations of relativistic axisymmetric MHD di sc winds h ave demonstrated that jets do 
indeed coUimate, with half openin g angles of 3-7° i Forth and Fendt .201Q,l . These authors 
set up a relativistic analogue of the lOuved and Pudritz (Il997h simulations, and emphasized 
flows far enough away from the central black hole that general relativistic effects could 
be ignored. The outflow has a mixed character — a combination of the Blandford-Payne, 
centrifugally driven wind with a toroidal pressure-dominated flow of the type discussed by 
iLvnden-Belll ( Il996h . Lorentz factors of the order 7^ = 6 are produced. The outflow has an 
interesting substructure, with a narrow relativistic jet towards the axis that is surrounded 
by a sub-relativistic flow launched from further out in the disc. This structure agrees with 
the observations of AGNs on the pc scales, as we discuss in Section 3. Finally we note 
that the toroidal fields that are generated in 3D simulations of jets seem to act as a sheath 
that protects the spine of the jet core from losing momentum to the surrounding medium 
llMignone et al.ll201Q) . 

The advent of 3D g eneral relativistic, MHD (GRMHD) codes (eg. iMcKiniievl j200d) : 
iMcKinnev etaP l2012ly ) has opened up the study of the role of black holes in the origin 
and evolution of AGN jets. As for the purely relativistic case, one of the primary ques- 
tions is to determine how relativistic the jets become, and to understand what maintains 
their stability and coUimation. It is generally thought that there are two reg imes for AGN 
jets: a magnetized disk surrounding a black hole, as described by the Bland ford and Pavn3 
jl982h picture, and/or the spinn ing black hole whose ergosph ere is threaded by magnetic 
field, as originally envisaged by iBlandford and Znaiekl (Il977h . Most theories and the sim- 
ulations suggest that it is the spinning black hole that produces the highly relativistic jets, 
while lower-speed outflows arise from the surrounding disk. Recent GRMHD simulations 
jMcKinnev and Blandford, 2009. ) model the initial threading magnetic field geometries for 
the spinning black hole as either a dipole or quadrupole configuration, while a thick torus is 
used to model an AGN disk. In addition to the need for significant spin of the central hole 
(the spin parameter of the hole is a/M > 0.4), their results demonstrate that the initial mag- 
netic configuration also plays a significant role. Jets produced by an initial dipole structure 
can achieve high Lorentz factors, {F > 10) and remain highly coUimated with an opening 
half-angle of 6j ~ 5° at 10^ gravitational radii. This is a stable jet with the substructure 
dominated by a stable m = 1 mode, probably excited in the underlying turbulent torus. 

However, this does not occur for a more complex, quadrupolar, initial magnetic geom- 
etry, where the jet does not reach steady realistic speeds {F < 3). The quadrupolar field 
leads to mass loading which considerably slows the jet, and makes it prone to disruption. 
This outflow has a weak and disorganized poloidal field and a much stronger (by factor of 
10-40) toroidal field. The authors speculate that these differences with the dipole jet may 
have something to do with the observed morphology of extragalactic jets; in particular, their 
classification into weakly relativistic "PR I-type" jets associated with AGNs in more clus- 
tered regions, and highly relativistic "FR Il-type" jets associated with isolated objects. The 
latter would be associated with the spinning black holes and the former with the underlying 
accretion disks. 

There are many similarities between these relativistic results and the non-relativistic 
work and theory described earlier. The stability of jets over a large variety of scales seems 
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assured for both regimes — they are able to run with stable substructure, dominated by an 
m = 1 mode, without falling apart. Both types of systems can remain very highly colli- 
mated, and this is related to their mass loading which is connected with their magnetic ge- 
ometry. Moreover, toroidal magnetic field dominates the non-relativistic jet structure, which 
again is a likely consequence of the theory we have outlined in this Section. 

In AGN systems, there is a concerted effort to understand jets on the smallest scales, 
which correspond to a few pc, in order to test theoretical models. The analysis described 
in this Section has shown that magnetic fields in jets are poloidal near the jet axis, and 
dominated by toroidal field farther out. This is a natural consequence of the acceleration 
and collimation mechanisms that should also apply to non-relativistic (eg. protostellar) and 
relativistic (AGN) accretion systems. We are now just beginning to be able to simulate the 
evolution of compact jets to larger scales (several pc) which allows direct comparison with 
observation. Simulations indicate that the field in protostellar jets becomes highly wrapped 
near the bow-shock of the jet (Staff et al. 2010), and this may have consequences for under- 
standing AGN jets as well. In the next Section, we turn to the parsec-scale structure of AGN 
jets, to connect observations with the theory and simulations presented above. 



3 Structure and dynamics of parsec-scale AGN jets 

3.1 Backdrop for current studies 

Radio studies of jets from Active Galactic Nuclei (AGNs) have primarily focused on parsec 
scales (with centimetre-wavelength Very Long Baseline Interferometry, VLBI, giving mil- 
liarcsecond resolution) and kiloparsec (kpc) scales (with arcsecond-resolution interferome- 
ters such as the Very Large Array and MERLIN), with relatively few observations probing 
intermediate scales. This means that, although the pc- and kpc-scale observations in some 
ways tell a coherent and self-consistent story, in others they seem to show discrepancies, 
with the connection between the behaviours observed on the two scales not always being 
clear. One reason for this may be that, although the parsec- and kpc-scale jets form different 
parts of a single structure, the observations may be determined by different contributions 
from "global" properties intrinsically related to the jet and "local" properties due to local- 
ized perturbations, interactions with the surrounding medium, variations in the magnetic 
field and density of the ambient medium, etc. 

Throughout most of the period in which AGN jets have been studied with VLBI, since 
the early 1980's, interpretations of observed features and behaviour have focused on the 
possibility that the observed jet properties are associated primarily with local agents, in par- 
ticular, relativistic shocks propagating in the jets. From a certain point of view, this was 
very natural, and was motivated by factors such as the visibly inhomogeneous appearance 
of the jets themselves (they are often dominated by distinct components) and the common 
occurrence of variability. A theoretical framework arose in which the jet components were 
primarily shocks propagating along the jet, with various processes related to these shocks 
giving rise to the observed very rapid variability. This picture seemed to receive support 
from the earliest VLBI polarization observations, which showed that an appreciable num- 
ber of AGN jet components have magnetic fields that are predominantly transverse to their 
jets — this was interpreted as reflecting the compression of an initially tangled magnetic 
field by a transverse shoc k, causing the field to become aligned in the plane of compres- 
sion dGabuzda et al .11 19921) . 
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More recently, various theoretical and observational studies have begun to explore the 
possibility that at least some of the polarization structures detected on parsec scales are 
associated with the intrinsic magnetic fields of the jets themselves. One area of such re- 
search involves theoretical simulations and observational searches for evidence of helical or 
toroidal magnetic fields carried by the jets, which are expected to form as aresult of the joint 
action of the rotation of the central black hole and its accretion disc and the jet outflow. The 
presence of toroidal or helical magnetic fields should give rise to characteristic transverse 
intensity and linear polarization structures across the jets, as well as transverse gradients in 
the observed Faraday rotation, due to the systematic change in the line-of-sight magnetic 
field across the jet. 

Thus, one of the challenges currently faced by researchers studying AGN jets is to try 
to discern which features are primarily due to the action of local agents, such as shocks or 
interactions with the ambient medium, and which are primarily associated with the intrinsic 
magnetic field and other properties of the jets themselves. This Section will review recent 
observational results for high-resolution (VLBI) observations of AGN jets, together with 
some related theoretical simulations, in the framework of this challenge. 



3.2 Estimates of VLBI-core B-field strengths 

The direction of the magnetic field B giving rise to synchrotron radiation can be deduced 
from the direction of the associated observed electric vector position angle (EVPA), if the 
optical depth regime is known: in the absence of Faraday rotation, B is perpendicular to 
the EVPA in optically thin regions and parallel to the EVPA in sufficiently optically thick 
regions (see, e.g., Gabuzda and Gomez 2001). Thus, it is relatively straightforward to get 
an idea of the overall magnetic-field configuration even from VLBI observations at a single 
frequency, as long as the effects of Faraday rotation are not too severe. However, it is not 
possible to estimate the strength of the synchrotron B field directly using observations at 
a single frequency: multi-frequency observations are required, and even then, somewhat 
indirect techniques must be applied. 

One of these is based on the frequency dependence of the position of the observed VLBI 
core in a Blandford-Konigl jet jBlandford and Konigj[l979l ). In this picture, the VLBI core 
is essentially like a photosphere, and corresponds to the surface where the optical depth 
T is equal to unity. The position of this r = 1 surface is located further down the jet at 
lower frequencies: r oc i^^^*^'', where r is distance from the jet base, u is frequency and kr 
is a parameter that, in general, depends on the spectral index and the laws for the decline 
of the electron density and magnetic field with distance from the central engine dKonigll 
Il98lh . In practice, the observed VLBI core includes emission from both the vicinity of this 
T = 1 surface and the innermost jet, but the essence of this basic theoretical picture holds. In 
contrast to the behaviour of the core region, optically thin regions — in other words, regions 
in the jet — are expected to coincide at different frequencies. By aligning optically thin 
regions observed simultaneously at different frequencies and deriving the relative positions 
of the observed VLBI cores, the resulting "core shifts" can be used to derive estimates of 
the core magnetic fields, if reas onable estimates for the Doppler factor, viewing angle and 
jet opening angle are available ( lLobanovlll998l) . This technique has only been applied to 
a handful of objects, but results obtained so far seem to give similar results for the cores 
of different AGNs, indicating that fcr ~ 1 in most cases (which corresponds to conditions 
not far from equi partition) and typical B fields a t a distance of 1 pc from the jet base of 
about 0.15 Gauss dO'SuUivan and GabuzdcJl2009ah . Although these estimates are somewhat 
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Fig. 4 Ex ample of a derived dependence o f _B-field strength on distance from the jet base, approximately 
B oc r-l jO'Sullivan and Gabuzdall2009j) . 



uncertain, due to uncertainties in the estimates of the Doppler factor and viewing angle, 
which cannot be derived directly from observations, the fact that similar core-region B fields 
are derived for various different AGNs suggests that these estimates are reasonably correct, 
at least to order of magnitude. These same studies have yielded evidence that the magnetic 
field falls off with distance from the jet base roughly as r^^ (Fig.|4]l. 



3.3 Basic observational properties of the polarization of AGN jets 

Although a wealth of polarization structure is observed among the variety of AGNs that 
have been studied with VLBl, it is possible to identify certain basic characteristics of the 
jet polarization. The most important of these is that, in a substantial majority of cases, the 
iet polarization is orient ed close to parallel or perpendicular to the local jet direction (e.g., 
Lister and Homanll2005l) . The effect of Faraday rotation is usually modest outside the core 
regions, leading to polarization rotations of no more than 10 — 15° at centimetre wave- 
lengths (see also Section 13.5b . For this reason, this tendency for the jet polarization to 
be either parallel or perpendicular to the local jet direction was evident ev en in the earli- 
est VLBI polarization observ ations carried out at a single wavelength (e.g.. lGabuzda et aD 
ll992l : ICawthome et al.lfl993h . On average, there is a tendency for AGNs with relatively lu- 
minous optical line emission ("quasars") to display higher apparent superluminal speeds and 
to most often have primarily longitudinal jet magnetic fields, while those having relatively 
weak optical line emission ("BL Lac objects") generally display lower apparent jet speeds 
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and most often hav e predominantly transverse jet magnetic fields (e.g., iGabuzdaet all 19941 : 
iBritzen et al.1l2008h . 

Jet regions with oblique, or "misaligned", polarization (i.e., neither close to parallel 
or perpendicular to the local jet direction) that cannot be ascribed to Faraday rotation are 
sometimes observed, but are comparatively rare. Jets with extensive regions of longitudinal 
or transverse B field are observed. Little is known about the nature of "inter-knot" polar- 
ized emission located between bright jet features; VSOP space- VLBI polarization obser- 
vations have revealed some cases in which longitudinal inter-knot polarization, implying 
transverse B fields associat ed with this underlying jet emission, are observed (iGabuzdal 
ll999l : |Pushkarev et al.ll2005h . 

Appreciable transverse polarization structure is observed in many jets, as becomes clear 
through a vi sual inspection of 15-GH z maps produced by the MOJAVE VLBA monitoring 
project (e.g.. lLister and Homanl200"5l) . In some cases, "spine-i-sheath" polarization structures 
are observed, with longitudinal polarization (orthogonal B field) near the cen tral axis of the 
Jet and orthogonal polarizati on (longitudinal B field) near the jet edges (e.g. lAttridge et al.l 
l l999l : |Pushkarev et al.ll2005l : Fig.|5li. In other cases, the predominant polarization is orthog- 
onal to the jet, and is offset appreciably from the central axis of the jet. 



3.4 Interpretation of observations 

Jet knots whose dominant polarization is longitudinal, implying a transverse B field, have 
frequently been interpreted as transverse shocks that have compressed an ini tially tangled 
magnetic field so that it has become ordered in the plane of compression ( lLaindll980l : 
iHughes etal.l 19891) . while the presence of orthogonal polarization (longitudinal B field) has 
been ascribed to the action of shear that has enhanced this component of the field in a layer 
surrounding the jet. Transverse polarization structures have similarly some times been in - 
terpreted as indicating the combined action of shocks and shear (e.g., Attridgeetal ][l99i). 
In the framework of the empirical connection between higher/lower apparent jet speeds 
and predominantly longitudinal/transverse jet magnetic fields, it has been suggested that 
transverse shocks may form more easily in jet outflows that have lower intrinsic speeds 
jGabuzdaetalll 19941 : iDuncan and Hughe sl l 19941) . 

Alternately, all of these polarization structures can b e understood as possi ble manifesta- 
tions of helical magnetic fields carried by the jets (e.g., iLvutikov et al]|2005l) : for example, 
jets with predominantly orthogonal or longitudinal B fields could carry helical fields with 
relatively large or small pitch angles. The tendency for AGNs displaying higher/lower ap- 
parent superluminal speeds to be more likely to have primarily longitudinal/transverse jet 
magnetic fields is interesting in this connection. One possibility is that this could reflect a 
relationship between the pitch angle of a helical jet B field and the physical speed of the jet 
outflow: the ratio of the speed of the jet outflow to the rotational speed of the central black 
hole and accretion disc could be higher for quasars, leading to smaller pitch angles for the 
helical B fields threading their jets and a dominance of the longitudinal component of the 
helical field, whereas the low outflow speeds of BL Lac objects lead to higher pitch angles 
for their helical fields and a dominance of the toroidal component of the helical field. 

The appearance of a "spine-Hsheath" transverse polarization structure or polarization, 
implying a longitudinal B field offset toward from the central axis of the jet, could likewise 
come about if the jet carries a helical magnetic field, with the azimuthal component dominat- 
ing near the central axis of the jet and the longitudinal component becoming dominant near 
the edges. The different viewpoints advocated by these two different types of interpretation 
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Fig. 5 Example of "spine+sheath" polarization structure across an the jet of the blazar 1055+018 
<Attridgeet anil999h . The object is at a redshift z = 0.888, so 1 mas equals 7.78 pc. The vectors show 
the direction of the S-field. 



essentially reflect the challenge referred to above, of distinguishing between characteristics 
associated with the action of local agents and those associated with intrinsic properties of 
the jets themselves. 

Considerable progress has also been made recently in using the results of theoretical 
computations and simulations of AGN jets to generate simulated VLBI intensity, spectral- 
index, polarization and Faraday-rotation distributions that can, i n principle, be compare d 
with observations. Examples inc l ude the studies carried out by [Z akamska et al. 



iGracia et al](l2009l),lMimica et al.l(l2009l).lBroderick and McKinnevlt201Q)„Porth and Fendii 



201C llPorth et alUlOl ih and lClausen-Brown et al.l(l201 ih . lMimica et alJ ( l2009h and lMimica and Alovl 
20ld . l2012 h have considered various observational signatures of the presence of shocks and 



shock -heate d gas in rela t ivistic jets, including jets candying helical fields. Broderick and McKinnevI 
( I2OIOI) and IPorth et all ( 1201 ih have generated theoretical Faraday-rotation images based 
on computations and simulations of AGN jets with a helical magnetic-field component. 
|Porth et al. (2011) also considered the behaviour of the polarization angle as a function of 
the wavelength A squared in the (partially) optically thick core region, demonstrating devia- 
tions from a relation away from a limited range of relatively short wavelengths; this effect 
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appears to have been observed by lO' Sullivan and Gabuzdal (l2009tih . One limitation here is 
the need to use some scheme to extrapolate results for small scales to the scales accessi- 
ble to observations with VLBI. However, such studies are beginning to enable much closer 
comparisons with observations than was possible previously, and are being used to identify 
a range of characteristic observational signatures of various types of jets and jet magnetic 
fields. It is clear that the ability to reliably measure transverse structure in intensity, spectral 
index, polarization and Faraday rotation across AGN jets will be key to probing the nature 
of the jet magnetic fields. 



3.5 Observational evidence for helical or toroidal magnetic fields 

As is discussed in Section 2 above, helical jet magnetic fields are predicted by a wide 
range of theoretical models and numerical simu l ations (e.g. Blandford and Pavn3 ll982l 
lOuyed and Pudritzll997l,lKrasnopolskv et alJl999i . lFendtl2006i) , making observational searches 
for evidence for helical magnetic fields associated with AGN jets of considerable interest. 
The most promising approach in this area is searching for transverse gradients i n the Fara- 
day ro tation measure distribution across AGN jets, as was first pointed out by IBlandford 
( Il993h . The observed Faraday rotation is proportional to the integral along the line of sight 
of UeB ■ dl, where ne and B are the electron density and magnetic field in the region of 
Faraday rotation and the length element dl points along the line of sight toward the observer. 
Therefore, the magnitude of the Faraday rotation depends on the line-of-sight component of 
the ambient magnetic field, and the sign of the Faraday rotation is determined by the direc- 
tion of this component of the field. Therefore, if a jet and its immediate vicinity are threaded 
by a helical magnetic field, this field will give rise to a systematic transverse gradient in the 
observed Faraday rotation, due to the systematic change in the line-of-sight component of 
the helical B field across the jet. The detection of such transverse Faraday-rotation gradients 
is observationally challenging, since it involves a joint analysis of multi-frequency radio im- 
ages, and the observing conditions are not optimal on either parsec or kpc scales: on parsec 
scales, the jets are very compact, so that they are sometimes only marginally resolved in 
regions of strong polarization, even when observed with VLBI, while, on kpc scales, the 
jets are better resolved, but the effect is weaker due to the increased distance from the base 
of the jet. To obtain increased sensitivity to Faraday rotation, one can observe at longer 
wavelengths, however, this comes at the expense of a corresponding loss of resolution. 

Nevertheless, firm detections of transverse Faraday rotation gradients have now been 
made for a number of AGNs. The first report of the detection of a transverse Faraday gra- 
dient across an A GN jet, interpreted as evidence for a helical jet B field, was made by 
lAsadaetal.l(l2002h . based on Very Long B aseline Array (VLB A) obs ervations of the nearby 



AGN 3C273; this was l ater confirmed by Zavala and Tavlo^ (l2005b . and the evolution of 
the gradient analyzed bv lAsada et all ( l2008bir In th e meantime, tra nsverse Faraday-rotation 
gradients in several more objects were presented bv lGabuzda et alJ ( i2004L) . There have since 
then been at least ten further papers in refereed journals reporting the presence of transverse 
Faraday-rotation gradients across the jets of more than a dozen AGNs on parsec scales, in- 
terpreted as evidence that these jets carry helical magnetic fields (Asada et al. 20 08b. 2008a|. 
20101 : [Gabuzda et al.l2008l: iGomez et al ll2008l;IO'Sumvan and Gabuzdall2009bl : iKharb et alj 



2009l :lMahmud et al.]l2009l : ICroke et al.H201(]l :lKronberg et al.ll201lh : one example is shown 



in Fig. [6l A number of other tentative cases have been presented in conference proceed- 
ings or identified in previously published Faraday- rotation images in the literature (e.g., 
IContopoulos et al.ll2009l ; rZavala and Tayloill2003l 120041 . 
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Fig. 6 Example of a transverse Faraday-rotation gradient detected across the jet of the blazar Mrk 501; the 
contours show the distribu tion of the radio intensity at 18cm, and the colour scale the distribution of the 
observed Faraday rotation jCroke et al.l2010ft . The redshift of this object is 0.034, so that 1 mas is 0.66 pc. 

This suggests that, when the observing conditions enable their detection, transverse 
Faraday-rotation gradients across AGN jets are not uncommon on parsec scales. This re- 
sult is potentially of crucial importance for our understanding of AGN jets, since it provides 
direct evidence that these jets carry helical magnetic fields, which, in turn, has obvious impli- 
cations for the jet-launching mechanism. In addition, the presence of an appreciable ordered 
toroidal field component would imply that the jets carry current. It is therefore important to 
try to ensure that the requirements for the reliable detection of transverse Faraday-rotation 
gradients are well understood and applied. 

It has sometimes been suggested that apparently transverse Faraday-rotation gradients 
could come about due to gradients in the ambient electron density, and therefore not have 
anything to do with the presence of a helical jet B field. An interesting point here is that, if 
the viewing angle of the jet in the jet rest frame is not too far from 90° (the jet is viewed 
roughly from the side in the jet rest frame), the observed Faraday rotation associated with 
the jet's helical field should change sign across the jet; if the jet viewing angle in the jet rest 
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frame is not close to 90°, the jet's helical field will still give rise to a transverse Faraday- 
rotation gradient, but there may not be a change in the sign of the Faraday rotation across the 
jet. Therefore, the clear detection of a change in the sign of the observed Faraday rotation 
across an AGN jet can easily be explained by the presence of a helical field, but not by 
elect ron-density gradients; such sign changes have been observed for a nu mber of AGN jets 
re.g.. lAsada et al.ll2008al : lGabuzda et all2004Ll2008l : lMahmud et al.ll2009h . 

Other observational problems that could potentially hinder the detection of transverse 
Faraday-rotation gradients associated with helical magnetic fields, or distort their interpre- 
tation, include the limited resolution available to centimeter-wave length VLBI observation s 
and optical depth effects in the region of the observed VLBI core /Tavlor and Zav alall2010t) . 
The presence of non-monotonic and/or complex transverse Farad ay-rotation gradients in the 
core re gion have been directly dem onstrated by the simulations of lBroderick and McKinnevI 
( l2010h and lPorth and Fendtl tOld) . These simulations indicate the difficulty of accurately 
deriving parameters of the jet and the helical field it is believed to carry based on observed 
transverse Faraday-rotation gradients. Note, however, that the direction of the simulated gra- 
dients is correct when they are monotonic; since most observational studies carried out thus 
far have focused on the detection of monotonic transverse gradients and their direction, the 
results of these studies are not greatly affected by the difficulties indicated by the simulation 
res ults. 

iTavlor and Zavalal (l2010h have proposed that a transverse Faraday-rotation gradient should 
span three "resolution elements" (taken to correspond to three beamwidths) in order for 
the gradient to be considered reliable; however, no theoretical basis was provided for this 
criterion. It is therefore of interest to develop more systematic and objective approaches 
to determining the resolution necessary to reliably distinguish transverse Faraday-rotation 
gradients, and results from a number of such studies have begun to appear. For example, 
iMurphv and Gabuzdal ( l2012h have calculated transverse intensity, polarization and Faraday- 
rotation distributions and subjected them to convolution with beams of various sizes. The 
initial results of these studies indicate that the detection of transverse polarization struc- 
ture associated with helical jet B fields is much more robust to convolution with beams 
comparable to the jet width than is transverse intensity structure, although, of course, the 
convolution somewhat distorts the intrinsic polarization structure. Similarly, convolution of 
transverse Faraday-rotation gradients with beams that are comparable to or even larger than 
the transverse size of the jet reduces the magnitude of the Fa raday-rotation gradient, but 
does not destroy the gradient completely. iHovatta et al.l ( 1201 2h have recently presented the 
results of four-frequency VLBI polarization observations of roughly 200 AGNs, together 
with Monte Carlo simulations modeling the polarization uncertainties and how they are 
manifest in Faraday-rotation maps. These Monte Carlo simulations suggest that transverse 
Faraday-rotation gradients spanning as little as 1.5 beamwidths can be reliable, provided 
that the difference between the Faraday rotations on either side of the jet exceeds 3a. 

Th ese results are also consistent with the general relativistic MHD simulations of lBroderick and McKinnc] 
( l2OI0h . which directly demonstrate the presence of transverse Faraday-rotation gradients 
across regions that are only marginally resolved: transverse structures in their simulated 
Faraday-rotation distributions with intrinsic sizes less than 0.05 mas remain clearly visible 
even when convolved with a 0.9-mas beam. All this suggests that the minimum resolution 
required to reliably detect a transverse Faraday-rotation gradient across an AGN jet is con - 
siderably less than the three-beamwidth criterion suggested by iTavlor and Zavalal (l20IOl) . 
provided that the Faraday rotation measures measured on opposite sides of the jet can re- 
liably be demonstrated to differ at more than the 3a level. Of course, this concerns only 
the detection of the presence and direction of a transverse Faraday-rotation gradient; reli- 
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able measurement of the intrinsic magnitude of such a gradient would require much higher 
resolution. 

Another issue that is currently being investigated is the question of how to accurately 
estimate the uncertainties at specific locations in VLBI images; this is non-trivial due to 
the complex computational procedures undertaken in the mapping process and the strong 
correlation between neighbouring pixels, due to convolution with the CLEAN beam. In the 
past, it has been usual to assign the uncertainty for a pixel in images of the Stokes parameters 
I, Q and U to be equal to the root-mean-square deviation about the mean calculated for 
a large region far f rom areas of source emiss ion, but this is certainly an underestimate. 
iHovatta et al. 1 ( I2OI2I) and lMahmud et alj j2012h use Monte Carlo simulations to address this 
issue, motivated by the need to derive accurate estimates of uncertainties for individual 
locations in Faraday-rotation images. 

Another interesting question in connection with the detection of transverse Faraday- 
rotation gradients across AGN jets is whether it is possible to observationally distinguish 
between a helical and a toroidal magnetic field. The point here is that the detection of the 
transverse gradient only demonstrates the presence of an ordered toroidal field component, 
without providing any information about the presence of a longitudinal field component. 
Both helical and toroidal jet B fields could give rise to "spine-i-sheath" transverse polariza- 
tion structures, and a rise in the degree of polarization toward the edges of the jet, as is com- 
monly observed. One means of distinguishing between these two possibilities is provided by 
the transverse intensity and polarization structure of the jet: a purely toroidal field can only 
produ ce symmetrical transverse intensity and polarization structures (e.g. Zakamska e t al.l 
l2008h . whereas a helical field can produce either symmetrical or asymmetrical transverse 
structures, dep ending on t he viewing angle of the jet and the pitch angle of the helical 
field (e.g., Murp hv et al EolO, 2012). Therefore, the joint detection of a transverse Faraday- 
rotation gradient and clearly asymmetrical transverse polarization structure would provide 
clear evidence for a helical, rather than a purely toroidal, jet B field. 

Finally, note that Faraday rotation associated with detected transverse gradients ob- 
served can be external (not occurring throughout the radiating volume of the source), if 
it is associated with regions of helical magnetic field in the outer layers of the jet, or in the 
immediate vicinity of the jet. 



3.6 Evidence for reversals in the direction of the toroidal field component 

One of the most unexpected results to come out of recent parsec-scale Faraday-rotation 
studies is the detection of reversals in the direction of transv erse Faraday-rotation gradient s 
across AGN jets , both with distance from the base of the iet lMahmu d and Gabuzdaji ; 
iHallahan andGa buzda (2008); Mah mud et all ^20121) and with time ( Mahmud et alJu 
This is at first glance difficult to understand, since the direction of the observed Faraday- 
rotation gradient is determined by the direction of the toroidal field component, which is, in 
turn, essentially determined by the direction of rotation of the central black hole and accre- 
tion disc, together with the direction of the net poloidal component of the initial "seed" field 
that is wound up. It is not physically plausible for the direction of the system's rotation to 
change on measureable time scales, and changing the direction of the jet poloidal component 
of the seed field would seem to require a change in the polarity of the field of the black h ole, 
which likewise, seems implausible, at last on measureable time scales. Bisnovatvi-K ogarJ 
( l2007h has suggested the possiblity that, under certain circumstances, torsional oscillations 
can develop in a jet, which could give rise to changes in the direction of the toroidal B-field 
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component; in this scenario, the direction of the observed transverse Faraday-rotation gradi- 
ents would reverse from time to time when the direction of the torsional oscillation reverses, 
and this rever sal pattern wou l d pres uma bly then propagate ou tward with the jet. 

However. lMahmud et alJ ( l200^ and lMahmud et alJ ( l2012l) have suggested that a simpler 
and more likely explanation is a magnetic-tower-type model, with poloidal magnetic flux 
and p oloidal current concentrated around the central axis (lLvnden-Bel]|l99d ; lNakamura et alj 
I2OO6'). Fundamental physics dictates that the magnetic-field lines must close; in this picture, 
the magnetic field forms meridional loops that are anchored in the inner and outer parts of 
the accretion disc, which become twisted due to the differential rotation of the disc. This 
should essentially give rise to an "inner" helical B field near the jet axis and an "outer" 
helical field somewhat further from the jet axis. These two regions of helical field will give 
rise to oppositely directed Faraday-rotation gradients, and the total observed gradient will 
be determined by which region of helical field dominates the observed Faraday rotation. The 
presence of a change in the direction of the observed transverse Faraday-rotation gradient 
with distance from the jet base could represent a transition from dominance of the inner to 
dominance of the outer helical B fields in the total observed Faraday rotation. The reversals 
of the Faraday-rotation gradients with time observed for 1803+784 ( iMahmud et alJl2009h 
could come about if the physical conditions in the jet changed such that there was a change 
in whether the inner or outer region of helical B field dominated the total observed Fara- 
day rotation. If this interpretation is correct, these observations of reversals of the direction 
of transverse Faraday-rotation gradients across AGN jets represent the first observational 
evidence that the jet magnetic field closes in the outer accretion disc. 

3.7 Evidence for the action of a cosmic battery 

Any transverse Faraday-rotation gradient can be described as being directe d either clockwise 
(CW) or counter-clockwise (CCW) on the sky, relative to the base of the iet. lContopoulos et alJ 
yOOQ) have reported a significant excess of CW transverse Faraday-rotation gradients for 
parsec-scale AGN jets, based on transverse Faraday-rotation gradients identified in maps 
from the literature. Considering the gradient closest to the VLBI core if two distinct regions 
with transverse gradients were present (as is described in the previous subsection) yielded 
29 transverse Faraday-rotation gradients on parsec scales , of which 22 were CW and only 
7 CCW, with the probability of this coming about by chance being less than 1%. This is an 
extremely counterintuitive result, since the direction of the helical field threading an AGN 
jet should essentially be determined by the direction of rotation of the central black hole and 
accretion disc, together with the direction of the poloidal "seed" field that is wound up, and 
our instincts tell us tha t both of these should be random. 

IContopoulos et al.1 (l2009h suggest that this seemingly bizarre result can be explained in 
a straightforward way via the action of a mechanism they call the "Poynting-Robertson cos- 
mic battery." The essence of this mechanism is the Poynting-Robertson drag experienced 
by charges in the accretion disc, which absorb energy emitted by the active nucleus and re- 
radiate this energy isotropically in their own rest frames. Because these charges are rotating 
with the accretion disc, this radiation will be beamed in the forward direction of their mo- 
tion, i.e., in the direction of the disc rotation. Due to conservation of momentum, the charges 
then feel a reaction force opposite to the direction of their motion; since the magnitude of 
this force exhibits an inverse dependence on the mass of the charge, this leads to a difference 
in the deceleration experienced by the protons and electrons in the disc. Since the electrons 
are decelerated more strongly, this leads to a net current in the disc, in the direction of ro- 
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Fig. 7 Illustration of how the action of the Poynting-Robeitson battery couples the direction of rotation of 
the accretion disc and the direction of the poloidal seed field that is wound up, leading to the creation of 
Faraday-rotation gradients asso ciated with the resulting helical B field that are CW on the sky, independent 
of the direction of disc rotation jMahmudll20 1 d) . 



tation. This current, in turn, gives rise to a net poloidal magnetic field whose direction is 
coupled to the direction of the current in the disc, i.e., to the direction of the disc rotation. 
This coupling of the disc rotation and the direction of the poloidal field that is "wound up" 
breaks the symmetry in the direction of the toroidal field component, and predicts that the 
observed Faraday-rotation gradients should be predominantly CW on the sky, independent 
of the direction of the disc rotation as seen by the observer (Fig. |7]l. In a "nested helical 
field" type picture such as that described in the previous subsection, the inner/outer regions 
of helical field should give rise to CW/CCW Faraday-rotation gradients; therefore, the ob- 
served excess of CW implies that the inner region of helical B field dominates on parsec 
scales. 

If this excess of CW Faraday-rotation gradients is confirmed by further studies, this 
will have cardinal implications for our understanding of AGN jets. The development of ap- 
proaches to carrying out objective searches for reliable transverse Faraday-rotation gradients 
for as many AGNs as possible is of interest in this connection. Indepe ndent of the physi- 
cal origin of the predominance of a particular orientation, the results of IContopoulos et aD 
12009) essentially rale out the possibility that most of the observed gradients are spurious, 
and not associated with the toroidal field components of the jets, since spurious (incorrectly 
identified) gradients should certainly be randomly oriented on the sky. If the Poynting- 
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Robertson battery does not operate sufficiently efficiently to provide the observed excess of 
CW gradients, then some other mechanism that can couple the direction of rotation of the 
accretion disc and the direction of the poloidal field that is wound up by differential rotation 
must be identified. 



3.8 Connection with kiloparsec scales 

Recently. IChristodoulou et al] ( 1201 2h (see also lGabuzda et al1l2012h searched the literature 
for transverse Faraday-rotation gradients across the kpc- scale jets and lobes of AGNs and ra- 
dio galaxies. This yielded six AGNs displaying continuous, monotonic RM gradients across 
their jets, oriented roughly orthogonal to the local jet direction. The most fundamental im- 
plication of this result is that transverse Faraday-rotation gradients that can most naturally 
be interpreted as evidence for toroidal or helical magnetic fields threading the jets are also 
present on kpc scales. This is an important result, because, if the jets of AGNs carry helical 
magnetic fields, these should be present on essentially all scales where the jets propagate, 
provided that the intrinsic field structure of the jet is not disrupted by interactions with the 
surrounding environment. ^^^^^^^^^^^^ 

Another intriguing aspect of the analysis of IChristodoulou et al] ( 1201 2h is that all six 
of the firm transverse kpc-scale Faraday-rotation gradients identified are oriented CCW on 
the sky, relative to their jet bases. This is on the verge of being a firm result statistically: 
based on a simple unweighted binomial probability function, the probability for six out of 
six gradients to be CCW by chance is about 1.5%. Howev er, the significance of this result 
is tenta tively supported by VLBA observation s of 3C 12 ("Couehlan et al.l l2010h . MrkSOl 
dCroke et al. 2010), and 1749-1-701 (Hallahan a nd Gabuz da 2008): in each case, transverse 
CCW Faraday-rotation gradients are found on scales of tens of pc. Further, all AGNs for 
which reversals of their transverse Faraday-rotation gradients have been ob served display 
CCW Faraday-rotation gradients further frorn the jet base (^ ahmud and Gabuzdal l2008l : 
iMahmud et al]l2012l : iReichstein and Gabuzdal I2OI0I : ICoughlan et al.l l2010). This growing 
evidence for a predominance of CW/CCW Faraday-rotation gradients on parsec/kpc scales 
can be interpreted as a consequence of a nested-helical-field configuration, as is described 
above, with the inner region of helical field dominating on parsec scales and the outer region 
of helical field dominating on kpc scales. This picture is consistent with theoretical studies of 
the jet launching mechanis m 1 Blandford and Pav"n3l 19821 : IContopoulos and Lovelacdl 19941 : 



IContopoulosll 1 995I : ISprui jl2010l) . which suggest that the field is effectively wound up only 
beyond the Alfven distance, which is ~10 times the radial extent of the outflow at its base. 
Thus, the theoretical models suggest that the magnetic field in the inner jet and the outer 
extended accretion disc will develop significant toroidal components on distances beyond 
-^10 AU and -^10 pc, respectively, from the AGN centres; this would predict a predomi- 
nance of CW Faraday-rotation gradients on observed scales less than ■^ lO pc and of CCW 



gradients o n scales greater than ~10pc, consistent with the results of IContopoulos et al 
2009) and IChristodoulou et all ( 120121) . Thus, although the results of IChristodoulou et al 



2012) must clearly be confirmed based on a larger number of AGNs, they have the poten- 



tial to considerably influence our understanding of the overall magnetic-field configurations 
of AGN jets. 

The analysis of IChristodoulou et al.l (l2012h suggests that the occurrence of reliably de- 
tectable transverse Faraday-rotation gradients is appreciably lower on kpc than on parsec 
scales: they were able to find reasonably clear transverse Faraday-rotation gradients in only 
six of about 85 objects. In fact, this is quite natural. There is an appreciable turbulent, in- 
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homogeneous component to the thermal ambient media surrounding the jets on kpc scales, 
which superposes a more or less random pattern over the systematic pattern due to the heli- 
cal fields (see the discussion in Section l43] l. This random component in the Faraday-rotation 
distribution may dominate in the majority of cases. Thus, it is natural that most of the ob- 
served Faraday-rotation distributions appear random and patchy, but the overall pattern due 
to the helical fields sometimes comes through. This suggests that, on average, it is easier 
to detect the systematic Faraday-rotation component due to helical jet magnetic fields on 
parsec scales, where the ordered inner field is more dominant. 

The additional complications imposed by the large-scale interaction with the jet and the 
ambient medium, and the structures (lobes, hotspots etc.) generated by this interaction, are 
discussed in the following Section. 

4 Interaction with surrounding medium and termination 

4.1 Introduction 

The best-studied regions of jet termination and the best examples of interactions of jets 
with the surrounding medium are furnished by radio-loud AGN, The parsec-scale AGN jets 
discussed in Section 3 generally feed into kpc-scale structures — kpc-scale jets, plumes, 
lobes and hotspots. The magnetic fields in these structures must exhibit some continuity 
with the fields set by the jet launching process (Section 2) and exhibited on the parsec scale 
(Section 3), particularly in the inner parts of kpc-scale jets. However, on the hundred-kpc 
scales of the jet termination and the downstream structures such as lobes and plumes, it 
seems likely that much of the initial structure has been erased. As is the case on parsec 
scales, the intrinsic polarization of synchrotron emission gives us information about the 
direction and degree of ordering of the magnetic field, though we always need to bear in 
mind that what we can measure is a projected, polarized-emission- weighted line-of-sight 
average of the intrinsic degree of polarization and the E- or B-vector angle. In addition, 
Faraday rotation both external and internal to the emitting medium has an effect on what 
we observe. A good deal is known about the polarization properties of radio-loud AGN on 
the largest scales, and we summarize the observational situation in Sections r4.3ll4.4l and l4.5l 
We then go on to discuss the relationship between observations and modelling in Section l4~6l 
and the situation in the comparable structures in stellar-scale outflows in Section l4~7] Before 
that, however, we discuss the methods by which magnetic field strength can be measured in 
AGN on these scales. 



4.2 Measuring magnetic field strengths on the kpc scale 

Unfortunately, synchrotron radiation on large scales can tell us little or nothing about the 
magnetic field strength, although this is clearly a key piece of information in understanding 
the dynamics and energetics of the jet and its environment, without significant (and possibly 
incorrect) additional assumptions. The methods discussed in Section [3^ can only be applied 
in regions with a well-known structure and a non-negligible optical depth, and kpc-scale 

structures are almost always optically thin at accessible frequencies^ 

The standard equipartition/minimum energy assumptions (Bu rbidg3ll956l) are widely 
used, but even at best only provide us with a plausible order of magnitude estimate of the 
magnetic field strength: there is no a priori reason to suppose that the plasma reaches the 
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minimum energy or equipartition conditions, and even if it does, these assumptions require 
us to guess the appropriate values for the fiUing factor and the fraction of energy in non- 
radiating particles, and k in the standard notation for equipartition: 

Finally, long-standing practice in the radio-galaxy community is to carry out the integral 
in eq. ( I13l l over a fixed range in observe d frequency, not in electron energy, although the 
latter seems more ph ysically reasonable dMvers and Spangled 1 19851 : iBmnetti et al.lll997l : 
iHardcastle et al .1 ^1 99 8) and the difference s between the se two sets of assumptions can lead 
to significantly different results (Beck and Krausell2005b . 

If equipartition and minimum energy are unsatisfactory, what else is available? Another 
widely used approach that gives at least some constraints is to consider the properties of a 
medium external to the radio source, which may be easier to measure. For example, thermal 
bremsstrahlung in the X-ray from the environments of large-scale jets and lobes can give 
us measurements of the temperature and density, and therefore pressure, of the external 
environment. On the reasonable assumption that the pressure in a radio-emitting component 
cannot be much less than that in the medium in which it is embedded, we therefore have 
a constraint on the total energy density in the source. Combining this with the synchrotron 
measurements, we can derive values for the magnetic field strength. However, these are not 
unique (for a given pressure greater than the minimum pressure there are two possible values 
of B to choose from) and the method is still dependent on assumptions about k and 4> as well 
as field geometry. 

By far the best method of measuring magnetic field strength in extended components 
of radio-loud AGN is the use of inverse-Compton emission. The volume emissivity from 
the inverse-Compton process depends essentially on the electron and photon number densi- 
ties as a function of energy. If the properties of the photon field are known, a detection of 
inverse-Compton emission gives us a constraint on the normalization of the electron energy 
spectrum. We can use this in combination with a measurement of synchrotron emissivity 
to estimate the B-field in a very robust way: assuming we know the geometry from radio 
observations, the only model dependence comes from our assumptions about the shape of 
the electron energy spectrum, since typically we will be observing energetically different 
electron populations using the two processes. For very simple spatial and spectral properties 
of the electron populations, analytical estim ates of B from inverse-Compton detections can 
be derived (e.g. iHarris and Grindlavl[l979h : for more complex situations, numerical codes 
must be used to integrate the relevant equations (e.g. IHardcastle etaDll998l : I Brunetti et al.l 
[2001; Hardcastle et al. 2002). 

Inverse-Compton em ission from large- s cale components of r adio-loud AGN was first 
discovered in the X-ray ('Harris et al.l ll994l: Ipeieelson et aDll995h and X-ray observations 
dominate the current applications of this technique of measuring B. A complication of this 
is that various other processes, including thermal bremsstrahlung from the hot phase of the 
IGM and synchrotron emission from high-energy electrons, can produce X-ray emission 
from the radio sources and their environments. In the following subsections we will discuss 
the available constraints on magnetic field strengths using this method, bearing in mind these 
observational limitations. 
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4.3 Magnetic fields in kpc-scale jets 



Kpc-scale jets are often strongly polarized, as discussed in Section !?^ In the most powerful 
jets, the apparent magnetic field direction after correction for Fa raday rotation is n ormally 
along the jet, often following bends in the jet very closely (e.g. iBridle et"aDll994) . In the 
less powerful jets of FRI radio galaxies, there is often a transition between a parallel field 
configuration in the inner few kpc to a perpendicular field i n the centre of the jet furt her out: 
parallel fields may be observed at the edge of the jet (e.g. bridle and Perle vl ll984l) . These 
FRI jets are the jets that have been studied in the greatest detail to date, because their prox- 
imity and relative brightness allows sensitive measurements to be made in regions of the 
jet that are transversely resolved on scales from sub-arcsec to tens of arcsec (hundreds of 
pc to a few kpc, at typical distances). In addition, both a jet and counterjet are routinely 
detected in FRIs, and the assumption that the jet and counterjet are intrinsically symmetrical 
is a viable one, allowing the effects of relativistic beaming and aberration to be separated 
from the intrinsic (rest-frame) behaviour of the source. Detailed modelling of these FRI jets 
(e.g. iLaing and Bridle|[2002l : Laing et al]|2006h shows that the observed polarization struc- 
tures can be modelled in terms of a magnetic field that is locally random but anisotropic. 
The anisotropic component of the field is mainly a combination of toroidal and longitudi- 
nal components, with the toroidal component dominating both at the edges of the jets and 
at large distances from the nucleus. Globally ordered helical models for the magnetic field 
(cf. Sectionl3.5b have been ruled out in some of the best-studied cases ( Canvin et al. 200^ : 
iLaing et alj|2006h : if the field is initially helical, it must evolve away from such a configura- 
tion by the kpc scale. 



Measurement of B-field strength s in jets is more difficult; althou gh resolved X-ray emis- 
sion from jets is common (see e.g. iHarris and Krawczvnskill200^ for a review), there are 
two fairly serious problems in using them to determine B. The first is that some of the 
X-ray emission seen from these structures, p articularly in low-power objects, is unequivo- 
cally known to be synchrotron in origin (e.g. iHardcastle et"al]|2001bl) . Where synchrotron 
emission is dominant, we only have a (usually uninteresting) lower limit on B in the jet. 
The second is that the j ets of the mo st powerful objects are known t o be at le ast mildly rel- 
ativistic on kpc scales (|Bridle et al.iri994 : Wardle and Aaron 19971 ; IHardcastle et al. 19*9^ ; 
iMuUin and Hardcastlell2009h and this introduces additional model dependence into any es- 
timates of B, since the observed properties of the jet now depend on the bulk Lorentz factor 
r and the angle to the line of sight 9. A good deal of interest was trigge red by the dis- 
covery that the X-ray emission from the powerful source PKS 0637-752 dSchwartz et all 
1^00), which could not be explained by a one-zone synchrotron model, could be explained 
as inverse-Compton scattering of the CMB with a magnetic field close to equipartition, us- 
ing a Lorentz factor (F ~ 1 0) and angle to the line of sight close to th ose implied by VLBI 
observations in the nucleus jTavecchio et al.ll2000l ; ICelotti et alj|200lh . which would imply 
no jet deceleration from pc to hundred-kpc scales. However, this certainly does not seem 
to be the case for all such powerful quasars ('Hardca sflelbooeh and later work on the broad- 
band SED of the knots in the best-studied example of this cla ss, 3C 273, suggests tha t the 
data are inconsistent with an inverse-Compton model anyway l l Jester et al ] |2002l , l2005h . As 
yet there is therefore no likelihood of being able to determine large-scale magnetic fields 
from inverse-Compton emission in the jets of the most powerful objects, though in (bright, 
nearby) lower-power objects we may hope for inverse-Compton emission at higher energies 
(see Section [53] l. 
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Fig. 8 Complex polarization structure in the lobes and hotspots of the FRII radio galaxy 3C 20 at 8 GHz. 
Right: the overall structure of the source. Top left: the E hotspot. Bottom left: the W hotspot. Polarization 
vectors have length representing fractional polarization and are plotted at 90° to the direction of the polar- 
ization i?-vector, thus giving a representation of the magnetic field direction. Contours show total intensity, 
increasing logarithmically by a factor of 2 at each step. The greyscales shown for the hotspots indicate polar- 
ized inte nsity to give an idea of the complex filamentary structure that is visible in polaiized emission. Data 
from iHardcastle et'aj] j 19971) . 



4.4 Jet termination 



The blight compact regions at the terminations of powerful jets, known as hotspots, are be- 
heved to be shocks at which the jet material interacts with plasma already present in the 
lobes. They are well studied in the radio be cause of their high surface brightness, and polar- 
ization observations at high resolution (e.g. ICariUi et al.|[l999t) show them to have complex 
polarization structures (e.g. Fig. |8), presumably because of the complex hydrodynamics 
in this part of the source. There is a general tendency for the inferred magnetic field an- 
ele to be perpendicular to the jet direction on entry to the hotspot (e.g. iLeahv et"al]| 19971 ; 
IHardcastle et al. 1997), possibly because of field compression at the shock. 

Because the hotspots are often the brightest compact features in sources in which they 
occur, they would be expected a priori to be good sources of so-called synchrotron self- 
Compton (SSC) emission, where synchrotron emission provides the seed photons for inverse- 
Compton scattering (Fig. |9]l. Observationa l constraints suggest that the bulk flow through 
hotspots is at best mildly relativistic (e.g. iMuUin et al.l l2008) and so beaming is less im- 
portant. Early observations of X-ray emission from hotspots were consistent in many cases 
with the X-rays bein g SSC and the magnetic fiel d strengths im plied being close to equipar- 
tition aHarris et al.i.l994ll200d IHardcastle et aT]|200 131120021) . However, at the same time, 
a number of hotspots showed X-ray emission that was much brighter than the expectation 
from SSC at equipartition and (in some cases) was spectrally inconsistent with being SSC 
at all; this X-ray emission seemed lik e ly to be partly or wholl y synchrotron in origin (e.g. 
[Harris et al.lll998l ; Iwilson et al.lboOlh . IHardcastle etd] ( |2004|) . based on the large amount 
of hotspot data collected in the first few years of the Chandra mission, showed that it was 
plausible that the difference was controlled by t he luminosity, and hence the magnetic field 
strength, of the hotspots (cf. lBrunetti et al .l2003h ; the most luminous hotspots have magnetic 
field strengths high enough that strong synchrotron losses inhibit the acceleration of elec- 
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Fi g. 9 Synchrotron self-Compton X-ray emission from the hotspots of Cygnus A (3C4()5), as first reported 
bv lHarr is et al. ( 199 ^. Lower pa nel: contours show the 4.5-GHz radio emission from the whole source, from 
the image of Carilli et alj 419911) . Upper panel shows the same contours overlaid on smoothed Chandra data 
in the 0.5-5.0 keV range from the dataset of Wilson et al. (2006). Note the close morphological agreement 
between radio and X-ray emission. The level of the inverse-Compton emission implies field strengths close 
to the equipartition value in the hotspots. 



trons to energies at which X-ray synchrotron emission would be detected. If this is the case, 
then it is plausible that all hotspots, not just the subs et for which SSC has been detec ted, have 
magn etic fields close to the equipartition value (Hard castle et a l. 2004; Kataoka and Stawm^ 
120051) . Remarkably, therefore, it seems that there is a rapidly-acting (since hotspots are short- 
lived), small-scale physical process that acts to amplify magnetic field strengths up to values 
within a factor of a few of their equipartition values. These field strengths may be up to 20 
nT (200 /iG) in the brighest sources. 



4.5 Lobes 

Beyond the end of the jet and the jet termination shock, if present, are the structures fed 
by the jet. The most common structure seen is large-scale lobes of synchrotron-emitting 
material extending back towards the centre of the host galaxy, although extended diffuse 
plumes are seen in a minority of low-power radio galaxies. 
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The lobes are often strongly polarized, with the typical apparent magnetic field direction 
being parallel to contours of constant total intensity surface brightness (e.g. Fig.[8ll. Because 
of the large spatial scales, strong polarization and large sizes of magnetic field structures, 
lobes are the best places for detailed studies of Faraday rotation intrinsic to the source or its 
local environment. The change in polarization position angle Ad produced by a foreground 
Faraday screen is given by Ad = RM X^, where RM is the rotation measure, 

RM = C / ne B ■ dl (14) 



Thus in principle, with multi-wavelength observations to give the required coverage, we 
can measure both RM and the intrinsic polarization angle and so constrain the electron 
density and field strength in the Faraday screen. In practice, there are a number of complica- 
tions (e.g. Laing 1984). Faraday-active material (thermal electrons) inside the synchrotron- 
emitting region will induce depolarization and departures from a law; in fact, the absence 
of this effect gives very strong, albeit model-dependent, lim its on the possible densities of 
thermal electrons inside the lobes (e.g. iDreher et al.lll987b . Even if the Faraday screen is 
external to the lobes, its characteristic size scale must be fully resolved by the radio obser- 
vations, otherwise beam depolarization results and the problem of deter mining Ue and B 
becom es much harder (though in some circumstances it is possible, see e.g. lHardcastle et al] 
|201(/). When Faraday rotation external to the radio galaxy is measured, the derived RM 
values are consistent with being due to the hot phase of the intracluster medium, and with 
constraints on the external density the field strengths may be estimated. Observed RM val- 
ues range from ±100 rad m ~^ in typical radio galaxies to ±4000 rad in Cygnus A, at 
the centre of a rich cluster dPreher et al.l ll987) and the estimated energy densities in mag- 
netic fields are then comparable to the thermal energy densities of the plasma. The structures 
observed in Faraday rotation may also be u sed to constrain th e power spectrum of magnetic 
field fluctuations in the external medium jLaing et al.ll2008l) . However, in the absence of 
internal Faraday rotation, for which there is no conclusive evidence in any object on kpc 
scales, these estimates of B tell us nothing about the field strength or its variations in the 
lobes themselves. 

The lobes are expected to be particularly good sources of inverse-Compton emission 
from scattering of diffuse photon populations such as the CMB and the extragalactic back- 
ground light (Fig. [Tol l. Because their geometry and synchrotron properties are also rela- 
tively well constrained, because it is certain that they do not have relativistic bulk motions, 
and because there is no evidence in general for significant ongoing particle acceleration, 
they represent the best cases for modelling and hence inverse-Compton observations give 
the cleanest constraints on magnetic field strengths. Earl y observations with R OSAT and 
ASCA allowed some field measurements to be made (e.g. Feigelson et al.lfl995h . but again 



it was Chandra and to s ome ext ent XMM-N ewton that gave us comparatively large samples 
( iHardcastle et al.ll200l ICroston et al.,.2004ll2005l : iKataoka and Stawarzll200"5.) . These pro- 
vide a very clear picture: both the detections and the non-detections of inverse-Compton 
emission are consistent with a picture in which the characteristic magnetic field in the lobes 
is close to, but in general slightly below, the equipartition value. The actual field strength 
value depends strongly on the lobe size and power, but might be of order 0. 1-few nT ( 1 to a 
few tens of /iG) in large, bright sources, with the energy density in electrons being perhaps 
an order of magnitude higher. Thus it seems that the mechanism that ensures (approximate) 
equipartition in the hotspot regions may also operate on the much larger, hundred-kpc scales 
of the lobes. 
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Fig. 10 X-ray emission from the nearby FRII radio galaxy Pictor A. The co lour scale shows Gau ssian- 
smoothed Chandra counts in the 0.5-5.0 k eV range, us ing the data described bv lMarshall et alj J2010I) : con- 
tours are fro m the 1 .4-GHz radio map of IPerlev et al.l (1997). Note the X-ray emission (most likely syn- 
chrotron, see lHardcastle and Crostoij2005l) from the jet and hotspot, the excellent match between the bound- 
aries of the X-ray and radio emission, and the relatively constant X-ray surface brightness of the lobes, in 
contrast to the strongly varying radio surface brightness (contours are logarithmic, increasing by a factor 2 at 
each step). 



We can go slightly further in the best-studied systems, since the fact that the synchrotron 
and inverse-Compton emission can be resolved with Chandra and XMM -Newton means tha t 
we are not restricted to considering one-zone models, as first realised by Isobe et alJ (l2002h . 
[Hardcastle and Croston (2005) showed that in fact the magnetic field strength has to vary 
by a factor of up to 2 throughout the lobe of the best-studied lobe IC source, Pictor A 
(Fig.llOb. in the sense that the field strength decreases with distance from the hotspot. More 
recently, Goodger et al. (in prep.) have used inverse-Compton data to argue that the strongly 
filamentary structure seen in synchrotron observations of many powerful radio sources must 
be a result of strong small-scale variations of the magnetic field strength (if it were the 
electron energy density that was varying, there would be a stronger correlation between the 
synchrotron and inverse-Compton emission than is observed). Although this detailed work is 
at the limits of what is presently achievable with X-ray observations, it may be the only way 
in which the detailed relationship between field and electrons in the lobes can be studied. 

The limitation of the lobe inverse-Compton technique is that it can only really be used 
to study sources in which the diffuse X-ray emission is not dominated by thermal brems- 
strahlung from the IGM; i.e. objects in relatively poor environments. Among other things, 
this means that with the exception of Cen A (see below. Section 15.3b there is no inverse- 
Compton estimate of the magnetic field strength in the lobes of any low-power (FRI) source. 
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although in verse-Compton limits have been used to constr ain the particle content in such 
objects (e.g. Ijetha et al.ll2008l : lHardcastle and Crostorill2010h . 



4.6 Modelling of large-scale fields 

Analytical modelling of the dynamics of the lobes of radio galaxies has since its beginning 
focus sed on a pure hydrodynamical approach (e.g. IScheueil 19741: Kaiser and Alexandej 



1997 ) and nu merical modelling has often followed suit (e.g. ISasson and Alexandeill2003F 



Krausell200 5l). To some extent, this approach may be (retrospectively) justifiable by the 



inverse-Compton observations discussed above (Section I4.5t . which tell us that magnetic 
fields cannot dominate the large-scale dynamics. Nevertheless, numerical modelling involv- 
ing magnetic fields is important for three reasons; firstly, even synthetic synchrotron map- 
ping requires some idea of the magnetic field strength distribution, while reproducing obser- 
vations of polarization and RM/depolarization structures (see above) is impossible without 
a model of the magnetic field both in the radio source and the environment; secondly, while 
the magnetic field may not dominate the dynamics of the large-scale structures, it may still 
have a dynamical role to play, for example in suppressing Kelvin-Helmholtz instabilities; 
and thirdly, such modelling in principle allows us to test assumptions about the transport of 
the magnetic field structures imposed at jet generation and observed in the parsec-scale jet 
(Sections 2 and 3). 

T he earliest simulations of lobe dynamics involving magnetic field structure dClarke et al.l 

ll989l ; lMatthews and Scheuejll990h involved passive transport of the field, but more recently 
it has been possible to numerically solve the full MHD equations for semi-realistic den- 
sity, jet speed and magnetization regimes. Some compromises have genera l ly been neces- 
sary to allo w such simulations. F or example, the work of iJones et al l ( 1199^ , iTregillis et al.l 
(1200 ih and lTreeiUis et"ZI (l2004l) considered electron transport and radiative losses, but op- 
erated in a uniform environment. MHD simulations involving realistic larg e-scale atmo- 
spheres, matched to the information derived from X-ray observations (e.g. lO'Neill et al.1 
l2005l ; iGaibler et al.ll200^ ; lO'NeiU and J Jl2010h tend to neglect electron radiative losses 
and also, in earlier work, tend to consider rather low jet density contrasts. I ncreased compu- 
tation al power has started to allow the merger of these two approaches (e.g'Mendy gral et al] 
l2012h . Key results from this modelling to date include the broad reproduction of the lobe 
dynamics expected from hydrodynamical analytic or numerical modelling; the suppres- 
sion of large-scale Kelvin-Helmholtz instabilities and of the entrainment of external gas 
seen in pure hydrodynamic modelling (Gaibler et al. 2009), and the emergence of com- 
plex small-scale str ucture in the magnetic field and hence the synchrotron emissivity (e.g. 
iTregillis et al .l200lh . 3D MHD simulations which can simultaneously capture realistic large- 
scale dynamics and have resolution good enough to be matched to the detailed observations 
of synchrotron and inverse-Compton emission in radio galaxy lobes, so as to help with 
questions such as the origin of filamentary structures seen in synchrotron emission (Sec- 
tion 14.5b can perhaps be expected in the next few years. Numerical modelling is also per- 
haps star ting to shed some lig ht on the origins of the field structures on large scales; for 
example, iGaibler et all ( 120091) put in as an input condition a helical field as observed in 
parsec-scale AGN jets (Section 3), and show that this initial field is amplified by shear pro- 
cesses to reach much larger v alues than expected from flux conservation. On the other hand, 
iHuarte-Espinosa etal ] (l201lh are able to reproduce many features of real radio galaxy lobes 
with a field which is initially random (in both the jet and the environment into which it 
propagates). It will be of great interest to see whether high-resolution simulations of large- 
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scale jet-environment interactions can give predictions which would allow us to distinguish 
between different models of the small-scale jet magnetic field structure. 



4.7 Large-scale magnetic fields in stellar-scale outflows 

The study of the terminations and the equivalents of lobes in stellar-scale outflows such 
as protostellar jets and microquasars is much less well developed. These structures require 
persistent jets, and so impulsive ejection events such as those seen in some microquasars 
are not expected to produce them. Where jets are persistent or at least recurrent on short 
timescales,, t here is evidence in a few well-studied microquasars for lobe-like synchrotron 
emission fe.elMirabel et al.ll 19921 : lFabrikall2004l : [Hardcastlell2doi) and even both lobes and 



hotspots llSoria et al.ll2010h . as well as indirect evide nce for the e xistence of lobes whose 
synchrotron emission is not seen in some other cases dOallo et al.| [2005). The lack of lobes 
in gene ral in these s ystems has been attributed to the low-density environments in which they 
reside ( lHeinzll2002h . Even when found, though, these lobes are generally faint, so that even 
polarization measurements are generally not possible, and techniques such as magnetic field 
strength measurements by inverse-Compton emission are completely out of the question, 
even setting aside the complex technical challenges of doing so in the presence of an often 
strong central X-ray source and of a complex and poorly constrained photon field. Authors 
wishing to estimate magnetic field strengths use the equipartition assumption, despite its 
limitations (Section l4~2l l. 

Turning to protostellar outflows, models for these systems fall into two general cate- 
gories: a jet-driven bow shock picture analoguous to the dynamics of FRII radio galaxies, 
and a wind-driven shell picture in which the molecular gas is driven by an underlying wide- 
angle wind componen t, such as is gi ven by the X-wind jCabrit et al.li 1997b . A survey of 
molecular outflows bv lLee et al.l ( l2000i) found that both mechanisms are needed in order to 
explain the full set of systems observed. The advent of Spitzer observations of protostellar 
jets has led to spectacular advances in our understanding of jet morphology, shock structure, 
and jet termination. For example, the jet in the classic T-Tauri star system HH46/47, imaged 
in the infrared with high resolution techniques (down to sub-arcsec levels), is bright and 
highly coUimated and has both a wide-angle outflow cavity as well as a coUimated jet that 
extend s through all scales fro m the protostellar scale up to the termination point on the bow 
shock JVelusamv et alj|2007h . showing a morphology which is remarkably similar to those 
of powerful radio galaxies. The jet is well coUimated with a width of 1-1.5" and length of 
10" (4500 AU). One of the most prominent features are bright "hotspots' that are coincident 
with the head of the jet and appear to be analogous to the hotspots of radio galaxies. Cool 
molecular gas is detected in the wide-angle flow that surrounds this jet. Similarly, observa- 
tions of the jet Cep E associated with a more massive forming star (with up to 4 Mq), using 
the same techniques, show shocked molecular hydrogen H2 emissio n — which requires th e 
existence of a magnetic precusor (C-shocks) in order to be excited dVelusamy et al.ll201 ih . 
Here again, the hottest emission comes from atomic/ionic gas produced at the hot spot at 
the bow shock of the jet. However, despite this important evidence for magnetization of the 
shocked region, and unlike the AGN case, there is as yet little direct evidence for the role of 
magnetic fields in the production of "hotspots" and lobes, since there are few observational 
diagnostics of field strength in this regime and synchrotron radiation is seldom obs erved. In 
the ca se of protostellar jets, the state of the art, as discussed above ( Carrasco-Gonza lez et al.l 
I2OIOI) is the detection of radio polarization, confirming a synchrotron origin for some of the 
detected radio continuum, but this is a detection of the jet, rather than of a lobe-like struc- 
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ture, and here again equipartition assumptions must be used to estimate a magnetic field 
strength. Next-generation radio and X-ray facilities will be required for the study of these 
systems to catch up with that of extragalactic large-scale lobes. 



5 New directions 

5.1 Magnetic field strengths and directions in protostellar jets 

The Expanded VLA (EVLA) and the extended Multi-Element Radio Linked Interferometer 
(e-MERLIN) will have an order of magnitude more sensitivity than their predecessors. Much 
larger parts of the radio spectrum will be measured as a consequence, and this will have a 
huge impact on our ability to measiu'e synchrotron radiation from protostellar systems. The 
one detected source HH80-81 is one of the most powerful TTS jets known. If the weaker 
jets turn out to have measurable synchrotron emission, this will herald a new era in the study 
of jets in general, and of their magnetic field structure in particular. 

A second major advance in understanding the launch of jets will come with the ALMA 
observatory. ALMA observations will allow one to resolve the 1-AU scale of nearby pro- 
tostellar discs for the first time. The important point is that ALMA will also be able to 
measure linearly polarized emission. It is expected that this will allow one to map the struc- 
ture of magnetic fields on the disc, and, in this way, to study the structure of the disc field 
and its connection to the launch mechanism. 

The bottom line is that the era of real magnetic field measurements in protostellar sys- 
tems is very near — and this will help revolutionize our knowledge of how jets are launched 
and coUimated. 



5.2 Parsec-scale magnetic fields in AGN: Physics of the launching, propagation and 
confinement of the jets 

Observations with the Very Long Baseline Array over the past decade have made the im- 
portance of high-resolution, multi-frequency, polarization-sensitive radio studies with VLBI 
abundantly clear. The availability of spectral information can help identify locations of par- 
ticle re-acceleration and low-frequency absorption in the jets, and the Faraday rotation dis- 
tribution can provide information about both the thermal gas in the vicinity of the AGN and 
the intrinsic magnetic fields of the jets themselves. The question of whether there is an ap- 
preciable helical or toroidal component is of crucial importance for our understanding of the 
launching and confinement of the jets, and will be a key aspect of future observational stud- 
ies. Multi-frequency polarization VLBI observations of more AGNs are very much needed 
for these reasons. 

Faraday-rotation studies on parsec scales with high resolution are of particular interest; 
however, work of this kind has been hindered by the fact that increasing resolution requires 
observing at higher frequencies, whereas increasing the sensitivity to Faraday rotation re- 
quires observing at lower frequencies. One promising way to overcome this limitation is 
through space VLBI observations at centimetre wavelengths. The currently flying RadioAs- 
tron mission is likely to be of only limited use for this purpose, since it has limited imaging 
capabiUty due to the elongated orbit of the space antenna. The Japanese-led VSOP2 mis- 
sion would have provided much better imaging capability, but, unfortunately, this mission 
has recently been cancelled. It will be important to ensiu^e that any other future space- VLBI 
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missions include frequencies low enough (8, 5 and possibly even 1.6 GHz) to make these 
power tools for high-resolution Faraday-rotation studies . On the other hand, polarization ob- 
servations at 3 mm-1.3 cm sampling subparsec scales could also prove valuable for future 
Faraday-rotation studies, since these will probe regions with high electron densities and high 
magnetic fields in the innermost jets, where very high Faraday rotations may be reached. 

Another area in which future VLBI observations can play an important role is in rel- 
atively low-frequency multi-frequency observations with global arrays such as the VLBA. 
Such observations can help link up the parsec scales studied this far with VLBI and the kpc 
scales that will be studied with high sensitivity using the EVLA and e-MERLIN. One such 
project is a recently completed set of four- wavelength 18-22 cm polarization VLBA obser- 
vations of the 135 AGNs in the main MOJAVE sample, which typically provide information 
about the jet structures out to scales of tens of parsec. Observations of additional AGNs at 
these and longer wavelengths would be of considerable interest. It is thus extremely impor- 
tant to develop techniques for accurate intensity, polarization and polarization position angle 
calibration of such long-wavelength VLBI data. 



5.3 Kiloparsec-scale AGN: new approaches to field direction and strength 

The EVLA, and other next-generation radio telescopes with observing bandwidths of or- 
der GHz rather than tens of MHz, will have a great impact on the study of Faraday rota- 
tion in the large-scale polarized components of radio galaxies. Rather than requiring many 
monochromatic multi-frequency observations, it will in many cases be possible to make a 
good measurement of RM and depolarization with a single broad-band observation. This 
should allow the results obtained from detail ed studies of bright objects in rich environments 
jPreher et al.lll987l : |Perlev and Carilli|[l996h to be generalized to more typical systems. The 
sensitivity and resolution provided by the EVLA and e-MERLIN will allow the detailed, 
quantitative study of polarization structures in jets (Section [4.3l l to be extended to more ob- 
jects, to smaller scales (connecting 100-pc and pc-scale jets) and possibly to the large-scale 
jets in FRII sources. At lower frequencies, instruments such as LOFAR with the capability 
to measure polarization sensitively at low frequencies may allow us to probe the RM due 
to very tenuous, large-scale gas seen in front of the polarized outer regions of giant radio 
galaxies, giving constraints on Ue and B on these scales that are currently unobtainable in 
any other way. 

Progress in inverse-Compton measurements of B will also have to come from new obser- 
vational developments. Currently one problem with the method is the extrapolation between 
the properties of the 7 ^ 1000 electrons responsible for inverse-Compton scattering CMB 
photons into the X-ray band and the 7 ^ 10* electrons responsible for radio synchrotron 
emission at GHz frequencies. This situation will be improved significantly with the advent of 
next-generation low-frequency radio telescopes such as LOFAR. We can also hope, with less 
certainty, for imaging missions working at harder X-ray energies, such as NuSTAR, which 
will not only probe somewhat higher-energy electrons but, by operating above the cutoff 
energies for thermal bremsstrahlung, will make it much easier to study inverse-Compton 
emission from sources in rich environments. Another ambiguity i n models is the correc t 
modelling of the very low-energy electron population: as shown bv lColafrancescol ( l2008l) . 
ALMA may in principle be able to detect the Sunyaev-Zel'dovich effect from radio galaxy 
lobes, which is a very strong probe of the total number of relativist ic electrons and hence 
of the low-energy cutoff (although see iHardcastle and Loonevll2008l for some caveats). Fi- 
nally, inverse-Compton studies at very high energies have some promise. Fermi recently 
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made the first detection of what appears to be inverse-Compton scattering of the CMB to 
~ GeV energies, in the nearby, and therefore very large in angular size, radio galaxy Cen 
A jAbdo et all201(t) : taken at face value the detection implies a magnetic field strength in 
the lobes that is again close to the equipartition value. Fermi can resolve only a very few 
nearby and (necessarily) atypical radio galaxies, but even limits from this process may be 
useful. At higher energies still, existing detections with TeV imaging instruments such as 
HESS and M AGIC place limits on inverse-Compton scattering (predominantly of starlight: 
IStawarz eTal. 2003) by the TeV electrons in the jets of low-power radio galaxies: for exam- 
ple, existing HESS detectio ns of Cen A, whose synchrotron properties are very well known 
(e.g. iHardcastle et al.ll2007l and refs therein) already require that the magnetic field cannot 
be too much lower than the e quipartition value, although detailed modelling is complex 
l lHardcastle and Crostonll201ll) . The resolution and sensitivity that will be achieved with in- 
struments such as the planned Cerenkov Telescope Array (CTA) will actually allow imaging 
of the TeV inverse-Compton jets in objects such as Cen A, if _B « Beq, and so will permit 
mapping of the B-field in an X-ray synchrotron jet for the first time. In the very long term. 
X-ray telescopes with very high sensitivity (e.g., the currently planned Athena mission) may 
allow us to measure inverse-Compton emission from the lobes of Galactic microquasars and 
synchrotron-emitting protostellar jets (Section [4.7t . 



6 Conclusion 



In conclusion, this review has covered a vast amount of observational, theoretical and com- 
putational phase space with the purpose of drawing out and analyzing the importance of 
magnetic fields for the origin and structure of astrophysical jets. The underlying engines 
for all of these systems involves rotating magnetized objects such as stars, black holes, and 
accretion discs. Relativistic and even general relativistic MHD leads to similar conclusions 
about jet dynamics to those deduced (and in some cases observed) in protostellar systems. 
While measurements of densities, temperatures, and velocities is straightforward in the latter 
type of system, magnetic measurements are very difficult. Exactly the opposite is the case 
for AGN jets. Encompassing these two complementary worlds lies the complete picture of 
astrophysical jets which is still in a very exciting stage of exploration. It will be impor- 
tant therefore to keep probing and extending our ability to measure magnetic field strengths 
and geometries in these diverse systems, since these provide insight into the origin of one 
nature's most marvellous creations. 
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